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ABSTKICT 


TWO  techniques  for  testing  ccanpact  heat  exchangers  are  refined. 

They  are  the  maximuE-siope  variation  of  the  single-blow  transient 
technique#  and  the  regulcir  periodic-  tcschnique. 

The  maximum-slope  method  is  improved  by  using  an  electronic  control 
to  produce  a  very  fast  heater  response,  thereby  generating  a  shcurper 
step  change  in  ciirstream  temperature.  The  design  method  for  the  elec¬ 
tronic  control  is  presented  along  with  an  experimental  verification  of 
its  function.  A  40  to  60  percent  reduction  in  the  heater  system  tire 
constcint  was  achieved.  This  produced  an  average  increase  in  exper:imental 
heat  transfer  results  of  4  percent,  for  four  giass-ceramic  matrix  sur¬ 
faces  tested. 

The  periodic  technique  is  improved  by  establishing  guidelines  for 

use  over  a  rcuige  0. 2  <  R  <  50.  Solutions  are  derived  to  include  the 

tu 

effects  of  wall  longitudinal  conduction,  both  finite  and  infinite.  The 
finite  solution  ia  incorpora,ted  into  a  data  reduction  program.  Experi- 
ntental  heat  transfer  results  obtained  using  the  periodic  technique  are 
presented  for  a  new  gl^s-ceramic  matrix  siurface  (Core  510}  in  the  test 
range  2.6  <  <  36.6.  ihis  range  wm  limited  by  test  rig  capabilif.ies 

with  respect  to  «dr  flow  rate  emd  healer  power,  not  the  test  techniqua. 
Periodic  technique  results  are  con^ared  with  meucimum-slope  method  results 
(3.6  <  <  21.3),  Tdxe  agreement  is  excellent. 

A  critical  cctnpariscn  of  the  periodic  and  maximum-slope  technique 
is  provided.  It  is  concluded  that  the  periodic  tecnnique  has  the  major 
advantages  of 

1.  A  very  wide  useful  range,  0.2  <  <  50. 

2.  Easily  generated  experimental  viaves. 

3.  The  accuracy  of  an  integral  technique. 

In  contrast  the  maximum-slope  method  possesses  the  advantages  of 

1.  Being  easy  and  fast  to  use. 

2.  Not  requiring  a  digital  computer  for  data  reduction. 
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1.  lOTRCffiUCTION 


•Jliers  is  currently  an  increasing  oi^asis  on  the  develo{®ent  of  the 
regenerative  gas  turbine  engine  for  vehicular  a^licatioi2.  One  of  the 
many  advantages  over  its  piston  engine  coi^titor  is  the  potential  for 
low  exhau^  gas  eoissions  (ssog) .  Several  regenerator  designs  incorpo¬ 
rate  a  periodic  flew  type  heat  exchemger  with  a  matrix  sxirface.  To  size 
the  matrix  surface,  designus  need  to  know  the  basic  flow  friction  and 
heat  transfer  charactei'istlcs  of  the  surface.  Heat  transfeir  character¬ 
istics  can  be  obtained  experimentally  from  an  aneaysis  of  the  transient 
response  of  a  test  core. 

«^ore  are  several  variations  of  fee  transient  technique  whife  in 
application  require  an  es^teL'isisiivtal  appeuratus  and  an  appropriate  mathe¬ 
matical  model.  The  term  math^atical  model  as  used  in  this  report 
consists  of  the  specified  idealized  system  including  fee  boundary  coitdi- 
tions  and  additionally  fee  differential  equations,  and  fee  soluclon. 

Heat  transfer  feeuracteristlcs  are  determined  by  matching  the  eiqperi- 
mentaliy  established  temperature  variation  of  fet  fluid  leaving  the 
matrix  with,  a  stailar  theoretical  result  trem  tiut  iafe^atical  model. 

The  mafeesatical  models  used  coesbnly  are  for  fee  st^  function  irput 
prtposed  by  lofee  [ij  ai^  others,  and  fee  periodic  function  input  brigi- 
eally  considered  by  Bell  and  Katz  [2j.  The%  models  require  feat  fee 
izraut  teeperatore  wav^ora  used  in  developing  the  theoretical  results, 
be  prodc^ired  by  the  eitPdiaei^tV'  This  is  efeecially  difficult  in  the 
case  of  the  step  functits). 

There  two  ways  of  correctizsg  for  errors  in  heat  transfer  results 
associated  with  non-st^  i  ^ts,  analytical  and  eaperincstal.  Analytical 
correction  involves  chai^iiq  tl»  test  ezroerisental  results  by  an  asount 
corresponding  to  the  deviation  fresa  the  true  step  function.  Kfelmayr 
[3}  has  dene  this  by  extending  Locke's  saxinas-sl^e  aefecd  to  arbitrary 
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upstze&m  fluid  te3|>erature  changes.  Kowever  his  findings  presented  in 

[3]  have  not  been  ei^ristentally  validated.  Eiqperiiaental  correctioa 
involves  ^veloping  a  test  rig  whidi  will  more  closely  approach  the 
step  function  input.  -This  approach  is  ttesirahie  since  a  soiaewhat 
j^culative  theoretical  "correction"  does  not  need  to  be  applied  to  the 
test  results. 

The  periodic  technique  was  considered  in  detail  by  Stang  and  Bush 

[4]  for  low  values  of  <i.e.  0.4  <  <  4.8).  The  technique  was 

found  to  be  oracticeil  and  further  developsient  to  extend  it  to  higher 

values  of  was  recoeaended. 

to 

^e  objectives  of  this  study  eure  to: 

1.  Design  and  build  a  heater  control  which  wien  connected  to  the 
Stanford  transient  test  rig  will  produce  a  step  in  t&spe-  ature 
which  Bore  closely  approaches  a  true  step  function. 

2.  £]q>erisaentaily  validate  the  cperation  of  the  new  heater  systea. 

3.  Bkaaise  possible  changes  in  heat  transfer  results  caused  by 
using  the  heater  control. 

4.  Further  develop  the  periodic  technique  so  that  it  w?ll  be 
usable  in  the  test  range  0.2  <  H.  ,  <  50. 

5i  Develop  a  sathesatical  aodel  for  the  periodic  techniqw  idiich 
will  acboont  for  the  influence  of  wall  loTsgitudina  1  ^  ccaduction. 

6;  Develx^  a  practical  data  reductios  procedare  for  the  periodic 
technique. 

2.  Cc^aze  h^t  transfer  results  cbtsii^d  using  the  periodic 
and  the  ^rigca-sicpe  aethod. 
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2.  GEKEP&L  Of  711E  TEST  RTO 


2.1  Introauctlon 

The  transiexxt  test  technique  uses  one  fluid  streaia  flowing  through 
a  test  aatrix.  The  general  schesatlc  diagram  of.  a  treinsient  test  rig 
such  as  that  used  at  Stanford  is  given  in  Fig.  2.1.  Details  of  this 
design  are  provided  by  Idiealer  [5].  Appendix  X  provides  detailed 
specifications  for  the  rig  components. 


Figure  2, 1  Transient  Test  Rig  sdxematic 


Fonctic^  of  tibe  -^^ious  el«sairts  axe: 

1.  Inlet  -  smooth  <nxt  vai^tions  in  fluid  velocity  and  t^^croture 
and  provide  a  saiboth  tinnsitioa  to  flow  in  a  duct. 

2.  Heaters  -  Provide,  an  ii^ut  te^eratore  response  for  the  air- 
streae  eiterii^  the  test  Mctixsi. 

3o  Heater  Control  -  Regulate  electric  pester  to  heaters. 

4,  Test  Section  -  Costa  ins  the  test  matrix. 
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5.  Flow  Ketor  ■»  Measure  accurately  the  air  flow  rate  through  the 
test  section* 

6.  Blower  ~  induces  air  flow  through  the  systen. 

7.  *  Pressure  and  tezi^rature  sensors. 

3.  Pitot  l^ibe  -  Two>dij»ensionai  velocity  traverse  downstream  of  the 
test  core. 

The  designs. of  the  inlet,  heaters,  test  section,  blower,  and  pressure 
and  taaperature  sensors  are  described  in  [5], 

When  the  heater  is  energized  it  functions  so  as  to  produce  either 
a  step^ch2Uige  or  a  cyclically  varying  teziperature  in  the  steadily  flow¬ 
ing  edxstream.  The  single-blow  transient  technique  employs  the  step- 
^ange.  The  periodic  technique  uses  the  cyclically  varying  tessperaturs. 

An  approximate  step-change  in  temperature  may  be  obtained  by  simply 
connecting  cr  disconnecting  an  electric  power  source.  This  tends  to 
produce  ah  exponential  type  of  temperature  change,  pointed  out  in  [5] 
and  shown  in  Secti<»i  3. 

Cyclic  teBq>erature  variation  may  be  obtained  by  usii^  a  DC  genera¬ 
tor  and  varying  the  ^ield  strength.;  Another  possibiUty  would  be  tP  use 
a  variable,  transfora^^  driven  lae^ahically.  These  could  rather 
unwieldy  and  viothK  arraixgement  using  electronic  circuitry  was  es^loyed. 
The  electranic  heater  control  circuit  is  described  in  the  following 
section. 


2.2  Heater  Control  Circuit 

The  circuit  in  Figs.  2. 2  and  2. 3  was  designed  to  provide  an  approxi- 
sate  step^^mge  in  tesq>erature  which  heis  a  short  rise-time  {<  0.05  sec). 
The  electronic  circ^try  of  Fig.  2.3  acc«»5>lishes  this  by  special  regu¬ 
lation  of  the  electrical  power  aipplied  to  the  heaters.  The  baisis  for 
1  he  design  wd  its  test  results  are  discussed  at  lei^th  in  Sections  3 
throv^h  6  and  Appendices  II,  IJI,  and  iv.  The  cirojdt  hM  the  addi¬ 
tional  capability  of  generating  a  cyclic  temperature  variation  when  an 
appropriate  c.  -lically-varyii^  resistor  is  used  (Periodic  liput.  Fig. 
2.3). 


A  list  of  the  coinponentR  used  in  the  heater  control  circuit  is 
given  in  Table  2.1.  Capacitors  and  are  banks  of  fixed  capacitors 
connected  using  manually  set  rotary  switches.  Schematics  of  the  switch¬ 
ing  diagrams  and  capacitor  banks  are  included  in  Appendix  I. 

The  descriptions  of  the  positions  of  switches  2  through  5  in  Fig, 
2.3  are; 

D  -  step-down  in  airstream  temperature 

U  -  step-up  in  airstream  temperature 

P  -  cyclic  variation  of  airstream  temperature 

N  -  normal  airstream  temperature  response  obtained  by  simply  turn¬ 
ing  ihe  electric  power  on-off. 

Nomenclature  for  the  silicon  controlled  rectifier  (SCR)  leads  is; 

K  cathode,  A  anode,  eind  G  gate. 

2.3  Operation  of  the  Heater  Control 

Wave  shapes,  showing  the  operation  ol  key  elements  in  the  heater 
control  circuit,  are  sketched  in  Fig.  2.4.  They  ray  be  used  as  an  aid 
in  understanding  the  operation  of  the  circuit,  or  as  a  guide  in  trouble¬ 
shooting  or  checking-out  a  newly  coiiStructed  circuit. 

The  basic  elect  ;ic  supply  is  117  volts  AC,  60  Hz.  This  sets  a 
limit  on  the  amo-int  <  *  electric  power  available  to  the  heaters,  which 
have  an  electrical  resistiince  of  5  ohms.  The  basic  AC  supply  is  full- 
wave  rectified  by  CRj^,  CR^,  CR^,  and  CR^  and  applied  to  the  firing 
circuit.  Fig.  2.3.  Tha  voltage  applied  to  the  firing  circuit  is  then 
"clamped"  or  held  almost  constant  by  zener  diode,  CR^.  The  effect  of 
this  clcimping  is  to  supply  the  firing  circuit  with  a  good  approximation 
of  a  DC  voltage  source. 

The  DC  voltage  is  applied  to  the  series  resistance-capacitance 
network  formed  by  R^  and  c^^.  The  voltage  across  C^^  will  be  a  simple 
exponential  function  of  time,  dependent  on  the  magnitude  of  the 
manually  variable  R2.  Thi.  voltage  formed  across  C^  is  applied  to  the 
emitter  of  unijvmction  transistor  (lead  1). 
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GE  C35B,  Silicon  Controlled  Rectifier 

GB  1H1694,  Rectifier  Diode 

GE  Z4XL22,  Zener  Diode,  22  volts 

GE  2N2646,  Unijunction  Transistor 

UTC  H51,  Pulse  Transformer 

330  n,  1  watt 

100  K,  Linear  Taper  Pot 

3300  a  5  watt 

470  0,  2  watt 

2000  0/  Linear  Taper  Pot 

39  X  n,  2  watt 

0.22  fjiF 

6  |JP 

0,5  jiP 
150  HP 
100  ^P 
80  {iP 
60  ^iP 

40  HP 
25  }4P 

20  |iP 
10  (iP 

SPST  Toggle  Switch 
Rotary  Switch,  4-position 
Rotary  Switch,  modified  to  function 
as  described  in  Appendix  I. 


EUECTRIC  POWER 
SUPPLY,  60  Hi 


VOLTAGE 
ACROSS  R3 


VOLTAGE  r  •  r 

ACROSS  CR5  L _ U. _ Ll 


i-inn 


VOLTAGE  y\  y\  y\  y 

ACROSS  C|  Lz.  lyLly  Lx  — 


GATE  CURRENT  I  I  I  I  I  I 
(6  SCR  TERMINAL)  _  R  L  i  l\  K 


VOLTAGE  ACROSS 
HEATERS 


Figure  2.4  Wave  shapes  of  the  Heater  Control  Circuit, 
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I5ie  unijunction  t3:ansistor  (OJT)  has  a  unique  electrical  charac¬ 
teristic  which  is  very  useful  in  this  application.  It  has  a  stable 
firing  voltage,  called  V^,  which  is  a  fixed  fraction  of  the  interbase 
voltage  V__  (lead  4  to  2)>  This  characteristic  is  shown  in  Fig.  2,5; 


Positive  Resistance- 
Characteristic 


Firing  Voltage 


Negative  Resistcuice 
Ch^kracteristic 


Emitter  Current 


Figure  2.5  Unijunction  Traujsistor  Emitter 
Characteristic  Cu^e 


As  the  emitter  voltage  V  (lead  1  to  2)  increases  from  zero,  the  emitter 
current  becooes  less  negative  and  then  slightly  positive.  Ihe  emitter 
in  this  region  behaves  much  like  a  high  passive  positive  resistor.  When 
V  reaches  the  firing  voltage  the  UJT  fires.  This  is  characterized 
by  tite  r  native  resistance  region  in  Fig.  2.5.  Firing  allows  the 
capacitor  C^  to  rapidly  discharge  through  one  winding  of  the  pulse 
transformer  Tj^.  A  pulse  voltage  is  then  induced  across  the  gates  of  the 
SCR’s.  Fig.  2,4  illustrates  the  gate  current  wave  shape.  The  resulting 
gate  current  triggers  the  forwwd  biased  silicon  controlled  rectifier 
(SCRj^  or  SCR^)  into  conduction. 

The  silicon  controlled  rectifier  (SCR)  is  a  unidirectional  semi¬ 
conductor  switch,  a  p-n-p-n  device.  In  its  blocking  state,  current  does 


not  flow  in  either  direction*  Hofwever,  if  a  saall  current  pulse  {typi- 

celly  3  volts  for  10  to  50  sticroseconds)  is  applied  froa  the  gate  G  to 

the  cathode  K  while  the  device  is  forward  biased,  tum-on  will  occur  and 

the  device  will  become  conducting.  Conduction  will  continue  until  the 

current  drops  below  about  100  bA  or  reverses.  In  the  reverse  biased 

direction,  the  SCR  acts  like  a  simple  dio^  and  blocks  current  flow. 

Because  of  the  very  fatst  tum-on  time  {caxi  be  considered  instantaneous), 

it  is  possible  to  trigger  an  SCR  into  conductii^  during  an  inteimsdiate 

point  on  the  AC  wave.  SCR  may  be  repeatedly  triggered,  cycle— after- 

cycle,  at  a  regulated  fraction  of  eacdi  AC  cycle.  The  result  is  a  control 

on  electric  power,  called  AC  lhase  Control.  The  circuit  shown  in  Figs. 

2.2  and  2.3  is  a  full-wave  ^ase  control  circuit  requiring  the  use  of  2 

SCR's.  Ref.  [6j  contains  a  broad  discussion  of  theory,  ratings,  and 

applications  of  the  SOI,  including  some  basic  AC  phase  control  circuits. 

The  unijunction  transistor  firing  voltage  is  a  f  inaction  of  the 

interbeise  voltage  V  .  The  specification  sheet  [7]  for  UJT  2N2646  lists 
SB 

this  function  as  typically 

Vp  =  0.65  Vgg  +  0.5  (2.1) 

The  interbase  characteristic  is  very  much  like  a  simple  resistance  if 

I_  ^  0.  The  firing  voltage  funchion  and  the  interbase  characteristic 

will  be  tised  in  the  attalysis  of  le.ter  sections.  By  regulating  V  ,  or 

BS 

the  time  required  for  the  voltage  across  to  reach  the  firing  tire 
or  phase  in  each  AC  wave  can  be  controlled.  The  resulting  AC  ^ase  con¬ 
trol  unit  will  be  used  as  the  heater  contml.  The  steady-state  output 
level  of  the  heater  coxrtrol  utilizes  the  time  required  for  voltage  across 
Cj^  to  reach  V^,  to  determine  its  magnitucte  and  is  controlled  mainly  by 
hn  initial  power  surge  (or  lach  of  it  in  the  step-down  condition), 
required  to  very  quickly  bring  the  heaters  to  steady-state,  uses  the 
variable  r«igoiaticn  of  v  .  This  later  tj-pe  of  elecrtrical  heater  power 
regulation  will  be  termed  the  AliSD  condition  in  the  sections  which  fol- 
Icn,  In  sui.'^iary,  the  heater  control  sharpens  the  step-change  by  applying 
em  electric  power  surge  (in  the  aided  step-up  case)  or  by  delaying  the 
lower  level  steady -state  power  (in  the  aided  step-dewn  case) . 
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3*  HEATER  ANALYSIS 


3.1  Description 

The  low  theraal  capacitance  airstream  heaters  are  described  fully 
in  [5].  A  brief  description  foiltws*  The  heater  system  consists  of 
four  banks  of  resistance  wires  electrically  in  peurallel  suspended  across 
the  airstream,  about  1.3  inches  apart  in  series  flow.  The  wire  is  0.003 
inch  diameter  alumel.  Each  bank  has  182  inches  of  wire  in  56  pcisses, 
each  3.25  inches  long.  The  wire  passes  ore  electrically  in  series. 

Thermal  expansion  of  the  wire  (about  0.01  inches  per  3.25  inch 
pass)  occurs  during  operation.  To  take  up  the  slack  each  pass  of  wire 
is  looped  over  an  extended  c2mtllever  spring  support.  Spring  deflection 
maintains  wire  tension  in  both  the  heated  and  unheated  states. 

Heater  wire  temperature  and  heat  flow  to  the  surrounding  airstream 
may  be  analyzed  using  the  thermal  circuit  in  Fig.  3.1. 


Figure  3.1  Thermal  circuit  for  Heater  Wires 


1 


With  respect  to  the  T„  variation,  the  air  temperature  Tac  can  be  treated 

H 

as  constant  amd  used  as  a  datum  (see  idealization  (4)  to  follow). 
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. . 


luedlizations  used  in  the  heater  analysis  are: 


1.  Wire  temperature  is  uniform  an'i  responds  as  a  single 

2.  fflieraal  resistance  is  a  function  if  air  flow  rate  only. 

3.  Itiemal  capacitance  C.  is  a  constant. 

H 

4.  llie  alrstream  teii^rature  change  in  lieating  20*F)  is  ssall 

cca^aured  to  T  {«*  3G0®F  above  the  ^dJ:  terr^rattare) . 

K 

Bie  thermal  resistance  for  convective  heat  transfer  from  the  heater  is 


defined: 


(3.1) 


3lie  thermal  'ipacltance  of  the  heaters  is  calculated  usii^ 


"  '“'h 


'3.2) 


Actual  values  for  and  are  shown  in  Fig.  3.2  and  are  calculated  in 
Appendix  II. 

Thermal  energy  rate  generated  by  the  electrical  resistance  of  the 
heaters,  is  plotted  in  Fig.  3.3.  The  figure  shows  this  energy 

rate  as  a  futKtion  of  the  voltage  Inpressed  by  the  heater  control.  la 
turn  the  voltage  can  be  regtilated  to  be  a  function  of  time. 

TSie  heater  in  Jig.  3. 1  has"  a  teaperatore  response,  T^,  ea^essed  in 
equation  form  as 


T  -T 
H  « 


=  q^(e) 


(3.3) 


The  energy  rate  (thermal  current)  flowing  through  the  convective  resistance 
is  deJi-vered  to  the  airstream.  Air  teiperature  rise  across  the  heater, 
ran  then  be  calculated  using  the  follcwing  energy  balance  on  the  air- 


stream: 


3.2  jcesposae  of  tiae  geaters  to  a  sisgle  Oa-Off  Voltage  Sas^y 

If  e  ccsstaiit  electric  volt^e,  cosisected  to  t3se  lieater  vires,  is 
siaply  svitched  oa  asd  off,  the  thezaal  esezgv  geaerated,  c  (8),  vili 

TO 

be  a  step  fui^rtirsu  She  solcticsi  to  £q.  (3.3),  ccsbiaed  vith  Eq.  {3.4), 
is 


Tj(0.ei  -  TjCO,®) 

r^{o,o)  -  TgCo,*)  " 


(3.5) 


»dsere  8^  ^  vhile  tbs  ^C9)  is  a  step  function,  the  afr- 

ze^csse  is  as  erposential.  Ihe  ^sastit^  8,.  is  called  the  tise 
coTistaat  cf  the  heaters  and  is  plotted  in  Slg.  3.2  for  the  losqied  case. 
As  Sjj  -»  O,  a  step-chance  ia  tessesatcre  is  approached. 

Jig.  3,2  preaente  ascther  tine— constant  carve  for  the  case  %ciere 
the  fear  basics  c£  heater  vires  axe  z»t  ccasidered  to  re^c^sd  as  a  lass. 
Ibe  beaters  are  cossi&red  as  roar  distinct  basics,  seqcgstially  pl«kced. 
in  the  airstrean.  £a^  bash  contsibates  abcat  23  percent  to  the  total 
tegperatore  rise  of  th%  air.  The  analysis  in  this  case  beetles  sore 
coe^Iicated.  Ihe  airstreass  tes^ratcre  7  is  a  tine^varying  feesrtien 
for  the  ^izee  dewnstreas  basics.  The  first  bank  of  heater  vise  is  still 
sid^ccted  to  caly  a  cocstabt  iZKcsing  airstrips  te^eratore.  Use  air- 
stzeaa  rejQ>SE&K  vrs  calculated  az^d  foosd  to  be  nearly  experiential. 

£g.  (3.3)  vas  fitted  to  the  calcclated  respoc^^  varying  8  until  a  goed 
fit  vas  c^Aained.  Ibe  values  obtained  for  the  effective  (tiK 
egatants)  are  plotted  in  Fig.  3.2  and  axe  seen  to  be  within  6  percent 
of  the  case.  For  this  reass.  it  is  believed  that  the  sisple 

li^ed  tseatnent  is  adc^jiate  for  design  purposes. 


Siespcase  of  Cie  Eeaters  ts  a  General  Tane-Varying  Voltage  Sc 


To  achieve  an  extrsely  short  rise  tine  in  the  heater  teegeraturg, 
the  sse  cf  a  tiae  varying  voltage  scrply  is  reqoired.  Regulatitgi  c£  the 
siqply  volt^^  in  tine  is  discussed  in  Secticsis  4  and  3  and  Appendices 
7TT  and  Z?.  Sie  thermal  energy  generated  in  the  heater  vires  will  also 
be  a  tiae-vasyiz^  fsnetixa^,  Sze  soluticsz  to  £c.  (3.3)  is  set  a  sisnle 
expo:»ntial  fus^tiea  and  rsqoires  the  use  of  a  ssnserical  aethod. 
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The  susierical  setzicd  was  the  ISM  Electrf^c  Circuit  Analysis 
Srogras  (ECAF).  Ihis  progrsa  was  written  to  be  used  hy  electrical 
eitgineers  in  the  design  and  analysis  of  circuits.  Detailed  infor&ation 
on  the  use  of  SCA?  say  be  found  in.  tl^  User's  Sfasoal  [s]. 

Use  of  £CA?  xequizes  ^ly  a  description  of  the  circuit,  a  list  of 
‘&e  elesents  in  each  branch  of  tbe  netsfoxk,  a  description  of  circuit 
excitaticm  oz  tise-vasying  elesents,  and  a  list  of  the  output  required. 

The  ti^~7arying  voltage  si^ly  to  tiie  heaters  can  be  converted 
into  a  tiae -dependent  current  source  for  use  in  £q.  (3.3)  by  using  Fig. 
3.2.  This  ccrrent  source  sust  be  ezttered  into  the  ECAF  prograa  in  the 
fors  of  actual  values  of  the  source  at  specified  tise  steps.  The  output 
of  the  SCAF  progret  was  »e  heater  wire  tesq^ature  in 

tom,  provided  tite  ri^t  ha:u3  side  of  £q.  (3*4)  idxich  was  then  used  to 
establi^  T^CO.O),  the  airstreaa  tesperature  downstreas  fron  the 
heater,  as  a  fosctisn  Of  tiae. 
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4.  STEP-UP  JCESIGN 


Step-up  in  temperature  occurs  when  voltage  is  applied  to  the  heaters 
and  the  airstrean  temperature  rises.  An  increase  of  about  20"F  is  used 
for  the  heater  control  design  and  in  the  esqperijBents.  If  the  heater 
voltage  applied  is  a  constant,  and  is  sii^ly  switched-on,  the  response 
c£  the  airstream  will  be  exponential  and  is  given  by  Eg.  (3,5).  Devia¬ 
tion  fr<xa  a  true  step-change  is  governed  by  the  magnitude  of  the  heater- 
wire  time  consteuit  0  .  As  0„  -♦  0,  a  step-change  in  air  teii5)erature 

is  approa<^ed.  The  value  of  0  has  been  partially  minimized  by  the  tise 

H 

of  small  diameter  wire  in  the  heater  £5].  Additional  decreases  in  the 
value  of  the  heater  system  effective  response  time  is  provided  by  the 
electrcsiic  control. 

If  instead  of  allying  a  constant  level  voltage,  corresponding  to 
a  steady-state  heating  energy  magnitude,  a  Icurger  voltage  is  applied. 

The  airstream,  of  course,  will  respond  faster  by  tending  toward  a  higher 
te^>eratnre,  process  of  applylr^  more  voltage,  than  is  required  for 
20®F,  is  termed  OVERSHOOTING.  Fig.  4. 1  shows  the  overrfiooting 
effect  cm  airstream  testpeiatuze.  The  part  of  the  step-up  response  curve 
of  interest,  is  the  initial  response  cheuracteristic.  A  ssnall  amount  of 
overshooting,  say  2X-Over8hoot,  will  reduce  by  about  60  percent  the  tlisg 
required  for  tlie  airstream  to  reach  63  percent  of  the  steady-state  value 
of  «  20®F. 

An  ideal  heate;t  control  would  initially  apply  a  volt^e  that  would 
produce  2X-Overshoot  or  more.  When  steady-state  is  reacdied,  this  applied 
voltage  would  revert  to  a  steady-state  level.  The  insulting  trai^erature 
wave  would  more  closely  resemble  a  true  step-change  in  airstream  tempera¬ 
ture,  with  a  SO  percent  response  achieved  in  only  25  percent  of  the  time 
required  by  the  normal  response  of  Fig.  4.1. 

TSie  heater  control  circuit,  presented  in  Figs.  2.2  cmd  2,3,  is 
designed  to  produce  an  overshoot  very  much  like  iiiat  inscribed  above. 

For  the  step-up  case,  SW2r  SW^,  SW^,  and  sre  set  so  position  o.  This 

has  the  effect  of  creating  initially  a  very  small  interbase  voltage 
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allowing  the  vni junction  transistor  to  fire  very  eeurly  in  each  AC 
wave.  The  net  effect  is  a  voltage  sut^e  to  the  heaters.  As  capacitor 
charges,  the  inteibase  voltage  increases  causing  the  voltage  applied 
to  the  heaters  to  taper  off  to  steady-state  level.  The  steady-state 
level  is  regulated  by  the  set  value  of  R^. 

An  cuialysls  has  been  done  to  determine  the  values  for  wd 
that  will  produce  the  desired  voltage  surge.  The  circuit-analysis 
technique  used  and  an  example  with  w  =  590  Ibm/hr  are  described  in 
Appendix  III.  The  result  of  the  analysis  at  each  flow  rate  is  a  tempera¬ 
ture  response  curve  for  the  heater  wires.  This  curve  is  very  nearly 
exponential  and  can  be  fitted  using  Eq.  (3.5),  varying  6  .  The  calcu- 
lated  time  constants  of  the  heater  wires  are  plotted  versus  air  flow 
rate  in  Pig.  4,2.  At  the  low  flow  rates  in  Fig.  4.2,  the  time  constants 
axe  less  th2m  at  the  higher  air  flow  rates.  This  is  the  result  of  more 
overshooting  ability  (lower  steady-'State  level)  at  the  lower  rates.  A 
comparison  of  Fig,  4.2  and  the  time  constant  curve  of  Fig,  3.2  indicates 
the  gains  realized  by  the  use  of  the  heater  control. 
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AIR  FLOW  RATE,  Ibm/hr 

Figure  4.2  Ai^d  Step-Up  Heater  Time  Conistants  from  Design  calculations. 
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5.  STEP-DOWN  ISSSIGN 


Step-down  in  temperature  occurs  when  the  voltage  applied  to  the 
heaters-  is  removed.  The  airstream  temperature  drops  and  the  response 
is  given  by  Eq.  (3.5).  Deviation  from  a  true  step-change  is  governed 
by  the  magnitude  of  the  heater-wire  time  constant  0  .  As  6  — ►  0, 
a  step-change  in  fluid  temperature  is  approached.  Hie  magnitude  of  6 

H 

has  been  partially  minimized  by  using  small  diameter  wire  in  the  heaters 
[5j.  Further  decreases  require  the  electronic  heater  control. 

Examining  the  step-down  response  curve  in  Fig.  5.1,  the  following 
is  noted.  The  response  of  the  airstream  initially  is  rather  quick  (steep 
slope)  until  about  75  percent  of  the  deflection  has  occurred.  If  the 
response  were  arrested  at  edsout  this  75  percent  point,  as  indicated  in 
Fig.  5.1,  a  response  curve  more  closely  resembling  a  step-change  would 
result.  This  could  be  accomplished  by  removing  the  voltage  to  the 
heaters,  and  shortly  thereafter  applying  a  lower  level  voltage  (enough 
to  produce  about  25  percent  of  ste^^Jy-state  heating). 

The  heater  control  circuit,  presented  in  Fig.  2.2  and  2.3,  is 
designed  to  remove  the  steady- state  volt2ige  and  a  short  time  later 
reapply  a  lower  voltage  to  the  heaters.  For  the  step-down  case  SW^^, 

SW^,  SW^,  and  SW^  eure  plcured  in  the  D  position.  This  has  two  main 
effects.  The  resistor  in  the  resistance-capacitance  network  which 
supplies  the  unijvmction  trcuisistor  emitter  voltage  is  increased  to  a 
maximum  level.  The  result  here  is  that  initially  the  UJT  fires  very 
late  in  the  AC  cycle  producing  almost  no  power  in  the  heaters.  The 
intexbase  voltage  then  decreases,  allowing  the  UJT  to  fire  progressively 
earlier  in  the  AC  wave  until  a  level  of  25  percent  of  steady-state  is 
reached.  Appendix  IV  includes  the  circuit-amalysis  technique  used  to 
calculate  the  response  of  the  heater  wires,  hence  the  airstream,  for  the 
step-down  case.  As  was  also  found  in  the  step-up  case  the  response  is 
nearly  exponential.  The  effective  time  constants  to  fit  Eq.  (3.5)  are 
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AIR  FLOW  RATE,  Ibm/hr 

Figure  5.2  Aided  Step-Down  Heater  Time  Constants  from  Design  Calculations. 
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6.  HEAT2R  COOTfiDL  PERFORMANCE 

6.1  Introduction 

Operation  or  the  heater  control  circuitry  was  investigated  experi- 
.entally,  as  the  resx»nse  predicted  by  the  design  in  Sections  4  and  5, 
needed  to  be  validated. 

With  no  core  in  the  test  section,  the  id.rstreaun  was  stepped  up  and 
down  (both  aided  wid  unaided) .  Temperature  response  o£  the  airstream 
was  monitored  using  thermocouples  and  a  fast  strip-cheirt  recorder.  With 
a  core  in  the  test  section,  chemges  in  the  heat  tramsfer  reznilts  due  to 
the  decreased  heater  time  constants  were  also  investigated.  Core  heat 
transfer  testing  tfiis  performed  using  both  the  elided  and  the  unaided 
3tep«ch^mges  in  fluid  temperature,  and  the  differences  in  results  were 
evaluated. 

6. 2  Heater-Response  Tests 

6.2.1  Unaided  -  The  first  set  of  experiments  performed  to  detezmine  the 
heater  response  had  no  core  in  the  airstream.  Since  the  thermocouples 
themselves  have  a  bead  diameter  which  is  l«u:ger  than  the  heater-wire 
diameter,  they  will  contribute  a  significant  leig  to  actual  airstream 
temperature  meeuiurement.  This  lag  effect  must  be  considered  because 
transient  airstream  temperature  is  being  measured.  R  thermal  circuit 
for  the  behavior  of  the  thermocouples  is  shown  in  Fig.  6.1: 


T^{0, 0) 


Figure  6.1  Thermal  Circuit  for  a  Thermocouple 
Suspended  in  the  Airstream 
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whera  is  sjiiere  diaroeter.  The  air  filjn  properties  were  taken  frora 
Bureau  of  Standards  Sources  [10]  and  are  evaluated  at  90»F. 

Ihe  heater- thermocouple  response  was  investigated  by  Wheeler  [5] 
shortly  after  the  construction  of  the  basic  test  rig.  Results  reported 
are  about  30  percent  lower  than  the  results  shown  in  Fig.  6.2.  The 
difference  is  likely  caused  by  the  methods  of  data  analysis  used.  Since 
the  thermocouple  output  record  is  not  a  simple  exponential  function,  see 
Eg.  (6.1),  the  results  of  a  curve  fit  will  depend  on  which  part  of  the 
experimental  trace  was  used  in  data  imalysis.  Results  shown  in  Fig. 

6.2  are  averages  of  several  points  read  over  the  middle  50  percent  of 
the  experimental  trace.  calculation  in  [5]  for  the  expected  heater 

wire  time  constant  (unedded)  is  about  20  percent  lower  than  the  results 
shown  in  Figs,  3,2  and  6.2,  Here  the  difference  is  due  to  different 
selections  of  both  gas  and  soLid  thermal  properties. 

To  determine  what  effect  the  therraocot^^le  time  constants,  jusc 
determined,  would  have  on  maximum-slope  heat  transfer  results,  the  ther¬ 
mal  circuit  of  the  thermocouple  in  Fig.  6,1  is  considered  with  the  T^ 
potential  source  input 

Tj(l,e)  =  T^(i,0)  4  [T^(l,«>)  -  Tj(l,0)]  (6,3) 


where  0  s  g  s  piau 

Ihis  "rac^*'  or  linear  input  is  used  because  it  resembles  the  eULrstream 
behavior  downstream  of  a  test  core,  especi2^.1y  at  the  time  when  the 
maximuto-slope  occurs.  For  the  thermal  circuit  in  Fig.  6.1 


T^(i,e)  - 


=  0 


(6.4) 


Solving  this  equation,  taking  the  derivative  dT^(G)/d9,  and 


comparing  with  dT,(l,  0)/^  yields 


d  •S  {S)/d8  — i— 

d  T^(l,9)/de  " 


1 

DG?tJ 


(6.?) 


I^e  above  ratio  can  be  interpreted  eis  the  ej^riisentally  established 

slope  of  the  air  downstream  of  the  test  core  (including  thenoocouple  lag) 

divided  by  the  actual  slope  of  the  downstream  air  (without  thermocouple 

lag).  Hie  maximum-slope  occurs  in  the  interval  0.3  KERU  <0  <0.5  DTRU, 

[14].  Additionally  frtm  meeisurements  presented  in  Fig.  6.2, 

9  <0.2  sec.  For  core  510,  DTAtl  >5.0  sec.,  so  the  exoonent  in  Eq. 

TC 

(6.5)  is  about  -10.  Consequently  it  is  concluded  that  thermocouple  lag 
has  a  negligible  influence  (much  less  than  0.1  percent)  on  maxiaisa-slope 
heat  transfer  results.  However,  leig  in  tka  heater  wires  does  have  a 
small  but  significemt  effect  on  the  heat  transfer  results  as  will  be 
shown  in  sections  following. 

6.2.2  Aided  -  After  the  thermocouple  time  constants  had  been  established, 
the  heater  systeo  qgerating  in  the  aided  mode  was  analysed.  Bctsed  on 
the  results  of  the  design  anal^is  in  Sections  4  and  5  and  Appendices 

III  and  IV,  the  heaters  can  in  the  aided  case  also  be  considered  to 
respond  as  a  single  exponential  function.  Ihis  allows  use  of  Eq.  (6.1) 
to  determine  the  reduced  heater  system  time  constant. 

Thermocouple  time  constant  0^  is  a  function  of  air  flow  rate  only, 
and  does  not  depend  on  the  heater  control.  The  vcQ.ues  plotted  in  Fig. 

6.2  for  the  unai^d  case  may  be  used  to  extract  tha  heater  tijae  constant 
in  the  aided  case.  Eq.  (6.1)  is  again  employed.  This  time  the  thermo¬ 
couple  trace  for  aided  step-changes  in  airstream  temperature,  and  0^ 
from  Fig.  6.2,  are  used  to  calculate  0  .  Resulting  heater  time  constants 
for  the  aided  case  are  plotted  in  Fig,  6.3  as  a  function  of  air  flow 
rate,  laie  design  prediction,  obtained  by  averaging  the  results  of  Figs. 

4.  2  and  5. 2,  is  also  included  in  Fig.  6. 3  for  coc5)arison. 
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Bie  results  sssarized  ia  Fig.  s.3  t2:st  a  ^seass  ia  tss 

heater  tise  constats  of  ^»at  60  pesceat  at  im  Ibs/hs  ts  40 

perceat  at  I0(X)  Ihss/hr  is  ps©&sc«a  hy  tm  beater  csaitrti  drcssitr^ 

Also  shown  in  tea  figt»  is  the  good  ^zee-Kst  hetace^  the 
cictioiaa  and  experinentallfr  detersdsed  reso'  ts. 

6.3  Sfrsotn  ca  Core  Heat  Srassfer  'Stests 

Heat  transfer  testing  was  perforaa^  co  fccr  test  SOSa,  SsSs, 

SOsC,  and  510.  sszis5>slcpe  sasfead  was  with  both  and 

unaided  step-ranges  la  fluid  tesnerature.  Fc^risestal  beat 
results  for  the  four  cores  are  plotted  la  Figs.  6.4  S.7,  5^5^ 

tlon  f^Ttor  results,  aim^  they  are  set  bv  ^  heat  traasSer 

teTciqi*  u^,  are  also  clotted  as  rtaiaed  for  the  aaa 

cases.  Bsis  re&s^aacy  icsares  feat  a?  errars  vers  in  the  air 

rate  calculations  by  prouiafag  a  cross  resde,  thros^  the  friction 
factor. 

Heat  trajisf cr  re^ilts  ehtaiced  using  the  alned  ste>-cbarge  were 
generally  higher  ran  correspessdiog  results  fsoc  tssts  with  the 
step-change.  Bteds  point  plotted  ia  Figs.  tsrsss^  6.7  is  ba^d 
eaperisentally  csi  tiie  average  of  two  ssd  tac  st^s-dbe=  is  air- 

stress  te^ezatore  «t*ring  the  cere.  AltheuT  r®  c=a:«e  is  resrOts 
can  be  disectly  re^  rros  fee  plots  faad  the  tabulated  results  in 

®^^ciaentail  variaticc  frees  rm  to  run  :2kes  it  di^icult  to  ses>- 
ing^iily  detezsine  the  saall  differed  as  a  fsjctioc  cf  air 
pr  Seysoids  nsser.  Artspta  were  ssae  to  test;  ersdit^sa  rusti¬ 
cal  in  both  ti»  ai^  and  ^r^dsd  cases,  ^suyi:^  rniy  &&  step- 

change.  Siis  was  accrr^liss^d  hy  perfesssiss  the  aids-?  and  rests 

bach-to-back  at  each  flcuj  rate. 

So  that  the  differences  in  heat  transfer  ressits  cu*sld  be  sassia^ 
sore  closely,  least— s^gjases  curse  fits  were  pw— — -.  -^3=  naasssss- 
sXoces  ana  the  ^  *s  of  each  core  ■versus  a— .»  K!^>i 

the  aided  and  unaided  ca^.  she  resslts  csf  iia  esr^  fits,  «?2lnafes 
at  ;^y:^lds  esdsers  in  Sie  respective  rashes  cf  festing.  are 
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CORE  SOSA 


in  Table  6.1.  From  the  table’s  listing,  the  average  change  in  heat 
transfer  results  for  the  four  cores,  when  using  the  aided  step-change, 
is  3.8  percent. 

Kohlmayr  [3, 11]  has  extended  the  naximum-slope  method  to  .include 
arbitr€ury  upstream  fluid  temperature  changes.  In  peorticular  he  reports 
tabulations  of  maximum  slope  MS  with  3  S  S  35  and  0.02  5  i  s  0.50, 
for  the  case  of  em  exponential  upstream  temperature  function,  changes 
in  maximum  slope  (/^),  predicted  by  [11],  are  plotted  in  Fig.  6.8  for 
5  s  s  35  aind  0. 02  s  ISO.  10.  The  pcurameter  I  is  defined  in  Table 
6.1  and  is  termed  a  deviation  from  the  true  step-change.  For  the 
values  of  I  emd  N  in  TcUDle  6.1,  [11]  predicts  that  the  change  in 

X*. 

maximum  slope  Z^IS  will  be  less  than  0.002,  a  change  under  0.4  percent. 
Using  the  sensitivity  curves  in  [12],  the  heat  tramsfer  results  will 
then  be  changed  by  less  than  one  percent.  However  the  experimental 
results  in  Table  6.1  indicate  an  average  change  of  3.8  percent  by  using 
the  aided  heaters  (which  do  not  produce  a  true  step-change,  only  a 
faster  response).  The  reeisons  for  the  disagreement  between  Kohlmayr’ s 
analytical  predictions.  Fig.  8.8,  and  the  experimental  results  reported 
here,  are  not  at  all  cleeir.  Further  work  on  this  question  is  needed. 
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HEAT  TRANSFER  DATA  CURVE  FITC,  FOR  AIDED  AND 
UNAIDED  STEP-CHAIKIE  CCMPARISCa? 


W  .  N  N  ^ 

Core  N_  f  MS  MS  MS  tu  tu  tu  ^  I 

K  «1 _ M _ %9.  '  '  ^  1  " 


"R 

Ibm/hr 

aided 

unaided 

aided 

unaided 

^tu'aided 

SOSA 

100 

107 

1. 138 

1.104 

0.034 

15.71 

14.73 

6.3 

.0087 

200 

214 

0.847 

0.334 

0.013 

8.04 

7.77 

3.4 

,0128 

300 

320 

0.713 

0.708 

0.005 

5.43 

5. 35 

1.5 

.0160 

400 

426 

0.631 

0.630 

0.001 

4.11 

4.10 

0.2 

.0168 

S0S3 

ICO 

107 

0.919 

0.896 

0. 023 

9.86 

9.33 

5.4 

.0131 

200 

215 

0.686 

0.678 

0.008 

4.89 

4.74 

3.1 

.0192 

300 

322 

0.578 

0.576 

0.002 

3.24 

3rl9 

1.5 

.0243 

SOSC 

100 

103 

0.699 

0.688 

0.011 

5.15 

4.96 

3.7 

.0265 

SIO 

100 

138 

1.160 

1.141 

0.019 

16.34 

15.66 

4.2 

.0081 

200 

275 

0.865 

0.850 

0.015 

8.53 

8.16 

4.3 

.0120 

300 

413 

0.728 

0.715 

0.013 

5.83 

5.58 

4.3 

.0151 

400 

550 

0.645 

0.632 

0,013 

4.^» 

4.25 

4.5 

.0175 

Pig-Jiria  5.8  Maxixti'jm-Slope  Correction  for  the  Exponential 
Temperature  Change  (from  Kohlmayr  [11]). 


7.  PERIODIC  lECHNIQUE 


7. 1  General 

The  regular  periodic  technique  for  determining  the  surface  average 
heat  transfer  coefficient  involves  imposing  a  cyclic  variation  in  tempera¬ 
ture  on  the  single  fluid  flowing  into  the  test  core  and  meeisuring  the 
temperature  response  of  the  fluid  flowing  out  of  the  core  after  steady- 
state  periodic  behavior  is  established.  The  measured  response  will  have 
an  amplitude  attenuation  and  a  phase  shift.  As  discussed  in  [4]  and  in 
more  detail  in  Section  7.3,  both  the  amplitude  attenuation  and  the  phase 
shift  of  the  fluid  temperature  response  may  be  used  to  obtain  accurate 
heat  transfer  results;  but  for  different  ranges;  thus  the  two  methods 
are  complementary. 

The  periodic  technique  has  been  Improved  by  the  follcwing  refine¬ 
ments  included  in  this  section: 

1.  A  solution  which  includes  the  effects  of  finite  wall  longi¬ 
tudinal  conduction  has  been  found  and  is  demonstrated  to  be 
’isable  in  data  reduction. 

2.  A  substantially  more  simple  solution  for  the  case  of  infinite 
wall  longitudinal  conduction  is  also  presented. 

3.  Guidelines  are  established  for  using  the  periodic  technique  tor 

0,2  <  <  50. 

tu 

4.  The  unceri;ainties  in  heat  transfer  results  ctssociated  with  core 
physical  prcpertxes  are  discussed. 

5.  Actual  test  results  for  a  new  test  core  (510)  are  presented 

for  2.6  <  <  36.6  and  are  co^>eured  with  the  aided  single- 

tu 

blow  results. 
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hU  lillltiui.  l>lll  >i)i» 


7.2  Solution  of  the  Governing  Equations,  Including  the 
Effects  of  Wall  Longitudinal  Conducticn 

7,2.1  Finite  Wall  Longitudinal  Conduction  -  The  governing  differential 
equations  for  the  case  of  finite  wall  longitudinal  conduction  are  derived 
in  Appendix  V  and  are 


*  * 

Tf  “  'i‘„  = 
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(7.1) 


(7.2) 


with  the  following  definitions 
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(7.3) 


The  wall  longitudinal  conduction  paraaeter  X,  and  N^^,in  conibination 
as  determine  the  magnitude  of  the  last  term  in  Eq.  (7.1).  The 

boundary  conditio  on  the  above  equati«i  is 

* 

T*(0,e‘>  =  sin  2t7-|-  =  sin-| 

with  t  =  ^ 

Zrr  Rc« 

The  solutions  to  Eqs.  {7,1)  and  <7.2}  are 
Fluid  Solution 


r*  =  «xp{-Z^Z)  I  sin 


(7.4) 
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wall  Solution 


*  * 

T*  =  exp(-fcj^Z)  Id-e^^)  sin(y  -  632 )  -  (7.5) 

The  terms  9Jid  are  functions  of  ^  and  ^^tu  must 

satisfy  the  following  set  of  equations 


2  /  3  2  2  2  \ 

*1  =  T  ■  -'vHi  *  'i  *  ^'1^2 '  '2 ) 


'2  *  ■  ^’‘tu(*2  ■  ^‘l'2  *  ^'l'2) 


These  constants  e,  and  e-  have  been  calculated  and  are  tabulated  in 
Tables  V-1  and  V-2.  For  data  reduction  the  Portrem  subroutine  C<snx:T 
in  Table  VII-3  was  vised  to  generate  the  tables  and  may  be  used  to 
calculate  e,  and  for  intermediate  values  of  ^  cind 

To  determine  the  surface  average  heat  transfer  coefficient,  Bq, 
(7.4)  is  used  and  the  temperature  response  is  meMured  downstream  of 
the  test  core  where  Z  =  From  Eq.  (7.4)  the  amplitude  attenuation 


D-  ^  ejqiC-e.N  ) 


(7.8) 


tJve  chase  shift  is 


Yf= 


(7.9) 


In  tile  special  case  of  X  =  0 


€,  - - r  aad  e, 

^  1  s- 


The  asciitude  attenuatior-  D_  and  the  psase  shift  are  ^o«n 

^  X  *  V 

graphically  in  Aiseadix  V,  Figs.  V-2  through  V-9,  as  functions  of 

.  f ,  and  X.  Sse  p^ase  shift  is  plotted  for  0  5  S  10, 

ist^4,  05X5  G.G2,  and  X  =  ••  The  acroiitwie  attenuation  is 

plotted  for  Q  5  5  50,  3  s  f  5  6,  O  5  x  s  0,02,  ass  X  =  ».  ©Je 

"tU 
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ranges  chosen  for  ^  cover  those  which  are  likely  to  be 

encotmtered  in  testing  using  the  two  methods. 

The  influence  of  X  on  experimental  amplitude  attenuation  is  shown 
in  Fig.  7.1  and  the  influence  on  experimental  phase  shift  is  shewn  in 
Fig.  7.2.  These  plots  indicate  when  the  effects  of  longitudinal  ccaiduc- 
tion  become  significant,  cos^red  to  other  experimental  uncertainties. 
They  also  can  be  used  to  estimate  the  magnitude  of  the  correction  that 
needs  to  be  applied  to  the  experimental  results  to  obtain  the  amplitude 
attenuation  and  phase  shift  that  do  not  contain  the  effects  of  longitu¬ 
dinal  conduction.  It  should  be  noted  that  these  plots  apply  only  when 
the  guidelines  for  testing,  set  in  Section  7.3,  aze  followed. 


7. 2. 2  Infinite  Wall  Longitudinal  Conduction  -  As  the  wall  longitudinal 

conduction  parameter  X  beccaaes  large  {X— ^  »),  the  tes^rature  of 

the  wall  T  becomes  a  functi<xi  of  t-itno  3  cmly.  Ihe  fluid  temperature 
w 


r-emadns  a  function  of  both  time  and  distance 


To  account  for  all  of 
the  thermal  energy  entering  and  leaving  the  wall,  the  entire  length  of 
the  heat  exchanger  must  be  considered.  Thermal  gradients  5r.  the  wall 
do  not  exist  and  this  is  the  only  means  of  applying  the  consertfatiwi  of 
energy  to  the  wall;  hence  the  resiilting  integral  Bq.  (7,10}  below.  Eq. 
(7,11)  for  the  fluid  is  also  presented  in  intergral  fora,  for  consistency 
reasons.  The  govemiig  equatiens  arc  derived  in  Appendix  V  ~nd  axB 
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The  solutions  to  ihe  abofVO  equatioafis  are 


~j-  =  crpC-CgZ)  sin  — 
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(7.12) 
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The  phase  shift  for  the  fluid  downstream  of  the  test  core  (Z  =  is 

plotted  as  the  X  ■=  <»  line  in  Figs.  V-2  through  V-5;  the  ao^litude 
attenuation  is  plotted  as  the  X  =  “  line  in  Figs.  V-6  through  V-9. 

Bell  and  Katz  [2]  in  1949  derived  an  exprossicii  for  the  amplitude 
attenuation  of  thsj  fluid  response  with  infinite  wall  longitudinal  conduc¬ 
tion,  for  S  2.  Considering  the  case  where  =  2  and  ^  =  1, 
the  calculated  amplitude  attenuation  is  0.4158.  Using  Eq.  (7.12)  with 
the  same  parameters,  produces  an  amplitude  attenuation  of  0.453.  Ihe 
disagreement  between  the  above  results  is  believed  to  be  caused  by  dif¬ 
fering  idealizations  used  in  the  .respective  analysis.  Ihe  analysis  in 
[2]  considers  the  effects  of  fluid  thermal  capacitance,  and  \ises  a 
mathematical  approxiisation  to  calculate  the  results.  The  analysis 
Included  here  in  Appendix  V,  considers  fluid  thermal  capacitance  to  be 
very  s^ll  (as  it  is  in  situations  involving  a  gas)  but  uses  no  mathe¬ 
matical  approximations.  It  is  believed  the  mathematical  approximation 
used  in  [2]  accounts  for  the  difference  in  the  results. 

7, 3  Guidelines  for  Using  the  Periodic  Technique 

Temperature  measurenent  uncert^d.nty  is  considered  to  be  the  largest 
source  of  error  in  the  periodic  technique.  Therefore  the  testing  should 
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be  done  in  a  !aa?5ner  that  vill  sinisdze  the  sensitivity  of  the  technique 
to  this  vincertainty.  A  tecperature  uncertainty  analysis  is  included  in 
Appendix  VI.  The  analysis  shows  that  the  sensitivity  of  the  periodic 
technique,  both  the  ar^Iitude  attenuation  and  the  {heise  shift  nethods, 
are  functions  of  and  For  a  given  there  exists  an  ^  vhich 

will  produce  the  best  test  results  (least  uncertainty  in  the  heat  trans¬ 
fer  results  for  an  uncertainty  in  tesqserature  neasurenent) .  If  an 
appropriate  value  of  ^  is  selected  emd  used  in  testing,  the  uncertainty 
in  heat  transfer  results  cam  be  held  at  about  3  pcxcent,  for  a  1  percent 
of  the  inlet  amplitude  uncertainty  in  teaperature  aeasurei^nt.  The 
values  of  rec<X3aended  for  testing,  and  the  method  tc  be  used  to 

reduce  the  data,  are  siicvn  in  Fig.  7.3.  This  figure  is  a  key  one  for 
use  in  setting-up  the  ea^risental  periods  of  oscillation.  Froa  the 
figure  and  the  analysis  of  Appendix  VI,  the  following  reccrssencations 
are  made: 

1.  <  1,  use  the  as^litude  attenuation  method  with  ^  <0.2, 

2.  is  S  7,  use  the  phase  shift  method  wxth  0.8  <  <  i,3, 

3.  N  S  7,  use  the  amplitude  attenuation  meteod  with  ^  >  2.9. 

* 

7.4  Sensitivity  of  the  Heat  Transfer  Results  to  Errors  in  Test 

Core  PhysiccOL  Properties 

The  periodic  transient  technique  requires  tliat  the  matrix  physical 
properties  such  as  density  p^,  specific  heat  c^,  and  then.tal  conduc¬ 
tivity  k  be  known,  in  the  periodic  technic[ue  if  is  Calculated  from 
w 

experimentally  obtained  amplitude  attenuation  and  phase  shift  Yf* 

From  an  uncertainty  analysis,  included  in  Appendix  VI,  and 
were  found  to  be  individually  sensitive  to  errors  in  But  because 

of  compensating  effects,  the  j — N  characteristic  is  not  affected  by 
the  errors  in  solid  density  (true  only  for  j  a  • 

^  p  ,  the  wall  density,  influences  the  Reynolds  Na  indirectly  through  its 
use  in  establishing  porosity  (flow  area)  gravimetrically.  see  Eg. 
(VI-9). 
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Figure  7.3  Rocojnmendod  iiondimonsional  Periods  Cor  Boot  Toot  Accuracy 


^  sssor  in  core  %peci£rc  heat:  does  affect  Ihe  heat  tratiofer  results, 
tut  net  the  nu^>er.  The  analysis  in  Appendix  VI  ^erws  that  the 

®®3^5to^®  o£  sensitivity  of  the  heat  transfer  results  to  errors  in  core 
spscific  heat  is  dependent  on  the  test  para=»ter  f,  The  sensitivitv  is 
plotiea  in  Fig.  VI-1-  The  flours  shenfs  that  if  the  periodic  teainique 
is  used  as  ruco'Wen&d  in  Section  7,3,  the  sensitivi^  for  the  anolitude 
attenuation  sathod  is 
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tu  tsi 

d  c  /c  “ 
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and  for  the  c^»se  ^lift  nethod  is 


^  5*-.  /K^ 
tu  tu 
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(7,,  16) 


(7.17) 


The  periodic  techni«aaa  then  is  very  se^'isitive  to  inaccuracy  in  tJie  solid 
specific  heat. 

Tin?,  periodic  te<duiique  is  about  as  sansitive  to  wall  longitudinal 

conduction  effects  as  the  saxisiua  slope  single-blow  techniqiKt.  Kail 

lor^itudinal  induction  effects  isost  be  considered  for  X  H  >  O.OL. 

tu 

Geiwrally  the  effect  cn  heat  transfer  results  is  less  than  10  percent. 
Therefore  an  uncertainty  in  solid  conductivity  of  ID  percent,  would 
produce  only  aboiit  a  one  percent  uncertainty  in  heat  transfer  resuits- 

7.5  Experiaental  Apparatus 

The  periodic  technique  uses  Ruch  of  the  saise  ecuipzsnt  as  the  single¬ 
blow  treuisient  technique.  The  wind  tunnel  test  rig,  originally  described 
by  Hheeler  [Sj,  has  been  nodified  as  described  laainly  in  Section  2  and 
Appendix  I.  For  fee  periodic  technique  fee  FUNCTI(»J  SWITCH  is  set  cn  ? 
{this  corresponds  to  the  P  positixxi  of  SW^,  SW^,  SW^,  and  SW_  in  Fig. 

2.3).  A  cyclically  varying  resistor  (stecbemical  drive)  is  connected 
to  the  heater  control  using  the  externail  reRIODIC  INPUT  plug.  This 
resistor  consists  of  two  resistors  in  series:  a  variable  resistor  used 
to  adjust  the  zero-level  (sdnintum  value)  of  fee  total  resistance  cuid  a 
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log-taoer  variable  resistor  «bich  is  driven  rv  the  can  ea  a  «>»li  vari- 
id>ie  speed  electric  »>tor.  !t$is  res^ilr  is  a  peri<^ic  variation  in  elec¬ 
trical  reaistante  «hich  is  traasfors»d  by  the  beater  control  into  a 
periodic  variation  of  the  ir^rat  of  eiectrical  ccwer  to  the  airstrees 
heater  wires.  Two  ssaJl  electric  cotofs  are  used  to  ©over  a  103:1  range 

of  period-of-escillaticjv  9  •  Cfee  aotor  has  a  torcus  of  52  in.  21s.  acd 

o  ~ 

a  speed  range  corresponding  to  7.5  <9  <  150  sec.  ihe  other  sator  has 

o 

a  torgce  of  4.0  in.  lbs.  and  a  spe^  range  corresponding  to 

1,5  <  S  <10  ssc, 
o 

The  resulting  airstreaa  tessezatsre  does  not  have  a  pure  sinusoidal 
rstperatnre  variation.  This  requires  the  use  of  a  Fourier  analysis  of 
the  te:5>erature  vaveforss  to  eatract  the  fundarasital.  To  illustrate  hew 
closely  the  innit  tesperattzre  wavefois;  resenbles  a  sissle  sinusoidal 
wave,  th“  test  run  for  Core  510,  n^pendix  VII,  with  %  =  226  is 
considered  here,  ©le  sr^Iitode  speetjnss  (i.e.  +  3^  )  ncosalised 

A 

to  the  first  haruffiiic  is  listed  in  Table  "/.  1: 

Table  7.1 

/donnse  siscmsi  tpe  mpur 
‘.'SSSiSSATOSE  «RVS,  SSSWaiTZSD 
to  ■nSE  FiaST  HaSHCSaC,  w  =  312  12ss/nr 

Magnitude 
1, 000 
0.053 
0,066 
0.028 

Clearly  there  is  a  fair  aiprox'nation  of  a  si^le  one-harsonic  sinusoidal 
te^rature  wave  inpuc. 

single  channel;  Honeywell  ElectroniJc  19,  lab  Hecorder,  is  used  to 
c^tedn  bo^  tha  upstxeas  and  dowxistreaa  tasperature  records.  To  do  this, 
the  recorder  is  ccamected  to  the  upstream  thersocoi^Xes  for  at  least  one 
period  9^,  and  then  quickly  switched  to  the  downstream  thensocoi^les  for 
at  least  one  period.  Using  the  svxtchijig  point  as  the  stcurting  point 


Harscnic 

1 

2 

3 

4 
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for  the  doKnstress  record  and  'backtarsckiog*  <s^  period  for  the  ns^szsext 

stsrtiag  poJxt^  the  tsrso  close- to-  siaaltaiatcss  tas«ratsre  wares,  relaxed 

is  the  data  seds^tica  progras  of  nzmeoilix  vU,  are  <g?t3i.aed. 

Ihe  thermocouple  s^sbsss  fas  upstreaa  ard  aosfsstxeas  each  cos- 

sisted  of  fear  series  ^esroco^les  csisg  a  sees  tenperfitase  zereny 

batii  as  the  cold  juartiss.  A  rocs  tenperatase  cold  jsr^tiicst  ves  used 

so  that  the  G-2  slllivolt  scale  cs  the  recorder  could  he  tssad  with  a 

saxissss  deflection  with  zero  offset  repaired,  rhe  tesraeratase  eapli— 

vas  held  at  shout  10®F  for  the  fluid  esterisg  the  test-  core. 

She  oeziod  of  oscillation  8  was  =sasu2rea  a  stcawat^h  and  a 
-  o 

count  of  the  saechanical  drive  «cles  on  the  variable  resistor. 

7.6  Results  of  Testing  Core  510 

7.6.1  Dessriptios  of  tSic  5est  Surface  -  Sse  surface  tested  was  Core 
So.  SIO  (Cemiag  Designatio  I.-»12).  It  is  a  slass-ccrasic  surface 
with  straight  triangular  pass^fcg.  She  core  rerves  a  daai  purpose. 

First  it  is  a  new  sarfars  of  the  ICKK)  cells/in^  class  ai^  prcfvides  basic 
heat  transfer  and  flow  friction  isfozsation.  Second  it  was  tested  using 
the  periodic  techaSgae  and  the  single-blow  technique  (ssaxiss  slope 
aethod)  ana  provides  a  direct  asErtPci^si  of  the  two  techniques.  The 
surface  geosetry  ano  coze  weight  are  described  in  Table  7.2: 


Table  7.2 


CeSE  510 


Si55Fi»CS  GSCeSSSuf  ASD  SSIGTT 


Cell  count,  H 
Porosity,  p 

Hydraulic  Diaaeter,  4rj^ 

Area  density,  a 

Cell  hei^V'-'iot^f 

h/4r. 

n 

Core  HasSff  K 


910  ceils/in'* 
0.716 

1.693  {io""  ft) 
1691  £t^/ft~ 
0.622 
134.5 
0.657  Iha 
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?.€.2  Sesalts  -  liis  hea&  trass£er  assi  £icv  frietica  test  resslrs 

are  plo*:t»a  ia  Fi^.  7.4  aad  are  listed  is  ^aSsle  VSZ-5.  2iie  £i^re 

deaszatraC'e?  ?exy  good  agreeaeat  betaeeen  tise  aided  sisgle>blov 

technisoe  (aaximars  slx^}  &e  pesioaio  fte  periodic 

tecsaioae  was  ssad  for  2-6  <  K.  <  36.6.  Shis  ranee  is  S.  wss  sot 
'  to  ta 

possible  for  tbe  TfaTiasaa  slooe  aeti^>d  bccasse  it  1*  limited  ^  S,  >  5. 

to 

Ste  effects  of  ^11  lasgit^^isal  caadactioa  oa  tbe  periodic  tachrir^ 
beat  transfer  zesslts  axe  ^sswn  oslv  for  tbe  hies  {low  ^  test 
poirts.  rest  of  ^e  test  points  were  virtsally^  anaffeeted  tbe 

wall  Icc^tadisal  ccodsctios.  frictios  f^tor  results  are  sot 

desesdest  oa  the  ty-pe  of  transient:  test  tectsiq:^  ssed  sisce  tastizg 
for  f  lswol#as  only  steady-state.  A  sssary  of  tse  basic  flow  frio- 
tior  a:^  heat  transfer  data  for  Core  510,  takes  froa  the  ssoctsed  csxves 
of  Fig,  7.4,  is  listed  is  2able  7.3. 

Heat  transfer  results,  are  abcot  percent  iewer  teas  'Q^oreti- 
cally  predicted  results  for  equilateral  triangular  tunes,  with  a  sosstast 
heat  flax  scundary  ocaditicn.  This  disagreesest  is  believed  to  be  catired 
primarily  by  a  ccabisatitgi  of  passage  zaaanlforaity  in  the  sastrix,  sos- 
triangular  shape  of  the  passag&s  (two-tsir^  of  the  come33  are  well 
rotizided},  and  ti^  heat  transfer  besEdary  tjosiditicsi  used  in  the  theory 
sot  being  well  dtmlicated  is  testing.  Friction  factor  results,  f,  are 
aboal  13  percent  lower  than  the  tlaory  for  equilateral  triangultr  passives, 
in  this  case  the  disagreesest  is  believed  to  be  caused  by  passage  a«»- 
ujiifossity  and  scn-triai^alar  shape;  frictassi  factor  is  not  affected  by 
the  heat  txsmsfer  boundary  c^mdition.  A  core  detailed  analysis  of  the 
above  effects  will  be  presented  in  a  future  report. 

Fig,  7.5  desosstrates  how  well  the  guidelinets  fzoQ  Section  7-3  were 
followed  in  tiie  Ore  510  tests.  The  experissstal  values  of  j  are 
plotted  for  the  corre^onding  values  of  figure  shows  that  the 

cental  values  of  the  nondisenslonai  period  ^  were  nil  within  the 
tolerance  lisd.ts  of  Fig.  7.3, 


Fiqiire  7.4  Core  510  surface  Characteristics  (ccasared  vith 
the  fully-developed  laainar  flow  ttieory  for 
equllateral-triacgular  passages) . 
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8.  COMPAniG  THE  PERIODIC  TECHNIQUE 
WITH  THE  MAXIMUM-SLOPE  METHOD 

8.1  Periodic  Technique 

The  periodic  technique  has  several  advantages; 

1 .  It  has  a  wide  useful  test  rang^,  0. 5  S  S  5C  (may  be  useful 
for  even  higher  values  of  . 

2.  It  is  conceptually  rather  siirple.  This  is  shown  in  the  analysis 
to  obtain  the  finite  wall  longitudinal  conduction  effects. 

3.  It  is  basically  an  integral  technique  vSnd  does  not  require 
extracting  derivatives  from  experimental  measurements. 

4.  Analytical  response  curves  are  not  required  since  the  heat 
transfer  results  are  computed  in  the  data  reduction. 

5.  Computing  time  is  short;  each  run  requires  about  1  second  on  an 
IBM  360/67. 

6.  Periodic  ten^ierature  waves  are  easier  to  generate  experimentally 
than  step  functions.  The  considerable  design  effort  included  in 
Sections  3,  4,  S,  and  6  could  be  eliminated  by  using  the  periodic 
twhnique^ 

Pi sadvantages  of  the  perio^c  technique  are; 

1.  It  requires  two  ten^erature  records. 

2.  A  digital  computer  is  required,  especially  for  the  Fourier 
analysis  emd  the  iterative  steps  used  in  accounting  for  wall 
longittidinal  conduction, 

3.  About  30  to  50  data  points  must  be  read  from  each  temperature 
record  and  this  information  prepared  for  tlie  computer,  requiring 
about  30  minutes  per  run, 

4.  A  suitable  period  of  oscillation,  0^,  must  be  detemdned  approxi¬ 
mately,  prior  to  testing. 


8. 2  Single?-Blow  Technique,  Maximum-Slope  Method 

Advantages  of  the  maximum-slope  method  include: 

1.  It  is  easy  and  fast  to  uses  The  experimental  temperature 
records  may  be  evaluated  and  prepared  for  the  computer  in  about 
15  minutes  per  run. 

2.  Computing  time  is  shorty  about  0.1  seconds  per  nui  on  an  IBM 
360/67.  Moreover,  the  computations  can  be  "hand  checked"  quite 
readily. 

3.  It  has  good  accuracy  for  5  <  <  50. 

Disadvantvges  of  the  maximum- slope  method  are: 

1.  It  requires  an  experimental  step-change  in  airstream  temperature 
which  is  difficult  to  achieve.  If  a  reasonably  good  step-change 
is  not  used  in  testing,  an  additional  correction  must  be  applied 
to  the  experimental  results  to  account  for  this. 

2.  Sensitivity  to  experimental  uncertainties  becomes  very  large 

in  the  0.5  S  S  5  range,  limiting  use  of  the  method  to 

N.  >  5. 
tu 

3.  Conceptually,  the  method  is  complicated.  Analytical  response 
curves  are  requited  in  the  data  reduction, 

4.  Experimental  derivatives  are  required, 

8.3  Similcirities 

The  periodic  technique  and  the  maximum  slope  method  are  very  similar 
in  the  following  ways: 

1.  only  one  fluid  is  used. 

2.  Matrij:  surfaces  can  be  tested. 

3.  No  wall  *-emperature  measurement  is  needed. 

4.  High  effectiveness  cores  can  be  tested. 
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5.  Wall  longitudinal  conduction  effects  on  test  results  are  about 
equal. 

6.  Sensitivities  to  uncertainty  in  core  specific  heat  are  about 

equal  for  10  <  <  50, 


54 


9.  CONCLUSIONS  AND  RECCMIENDATIONS 


9. 1  Conclusions 

1.  Kie  heater  control,  designed  and  constructed  as  a  part  of  this 
project,  functions  satisfactorily.  The  time  constant  of  the 
heater  system  is  reduced  by  from  60  percent  at  100  Ibm/hr  to 

40  percent  at  1000  Ibm/hr.  Some  of  the  otlier  desirable  features 
of  the  heater  control  are;  solid-state  device  with  no  mainte- 
nemce  required,  small  size  emd  economical,  control  by  a  small 
resistor,  flexibility  in  application  to  systems  requiring  more 
power  by  just  changing  the  silicon  control  rectifiers,  and 
cyclic  temperature  control  possibility. 

2.  Using  the  fast  response  heater  system  (aided)  caused  the  heat 
transfer  results  to  be  higher  by  an  average  of  3.8  percent,  for 
the  four  glass-ceramic  cores  tested. 

3.  Guidelines  for  using  the  periodic  technique  were  extended  so  it 

new  can  be  used  in  the  range  0.2  <  <  50.  Use  at  >  50 

tu  tu 

was  not  considered,  but  should  pres'-.nt  no  special  problems. 

4.  The  periodic  technique  was  generalized  to  account  and  correct 
for  the  effects  of  wall  longitudinal  conduction.  The  magnitude 
of  the  corrections  was  fotmd  to  be  about  as  large  as  for  the 
maximum-slope  method. 

5.  Hie  periodic  technique  was  used  for  2.6  <  N  <  36.6.  Very 

X-M 

good  agreement  was  demonstrated  with  the  maximum-slope  results 
obtained  using  the  fast  response  heater  system. 

6.  As  demonstrated  by  the  above  agreement,  both  the  periodic  and 
the  maximum-slope  techniques  give  good  results  when  properly 
applied. 
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9*2  Reconmendations 


i*  A  fast  response  heater  system  is  required  in  the  maximum  slope 
method.  A  heater  control  circuit,  such  as  presented  in  this 
report,  is  recatsmended  as  a  means  to  achieve  the  fast  response. 

2.  Ihe  periodic  technique  is  reconanended  ar  the  overall  best 
transient  technique  to  use.  A  wide  N^^-test  remge  as  well  as 
an  easily  generated  temperature  wave  are  two  major  adv£uit3ges. 

3.  To  resolve  the  differences  between  the  change  in  heat  transfer 
results  obtained  experimentally  in  this  project  due  to  a 
sheurper  input  step-change,  and  the  analytical  predictions  of 
[3, 11],  the  following  are  reconanended* 

a.  Perform  maximum-slope  tests,  intentionally  using  em 
extremely  slow-response  heater  system.  This  would  allow 
easy  and  accurate  measurements  of  the  input  temperature 
function  and  would  produce  laurge  effects  on  the  heat 
transfer  results. 

b.  use  a  geranetrically  uniform  test  core  to  eliminate 
possible  effects  or  passage  nonunifoxmlty  on  results. 

4.  the  effect  of  a  non- simple-sinusoidal  input  wave  cn  heat  trans¬ 
fer  results  using  the  periodic  technique  should  be  determined. 
One  way  to  experimentally  obtedn  a  more  reguleu:  sinusoidal  wave 
is  by  placing  a  low-N^^  matrix  in  the  airstream  between  the 
heaters  and  the  test  core.  !li1ie  matrix  would  serve  to  damp  out 
higher  harmonics. 

5.  The  sensitivity  of  the  periodic  technique  to  nonuniform  flow 
distribution  in.  the  test  matrix  should  be  analyzed. 
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APPENDIX  I 


GENERAL  TEST  RIG  INFOK'iATlON 

SpeGifications  of  the  Tegt  Rig  (see  also  Fig.  2.1) 

Table  II-l 

STANFORD  TRANSIENT  TEST  RIG  SPECIFICATIONS 

Inlet;  Belb&outh  with  a  pack  of  screens  and  an  egg-crate  flow 
straightener.  A  fan  is  used  to  mix  the  air  before 
entering  the  inlet. 

Heaters:  Four,  banks  of  0.003  in.  dia.  alumel  wire  with 

a.  Maximum  heat  generation  rate  -  2.37  Btu/sec 
(=:  2500  watts) 

b.  Electrical  resistance  -  5  ohms 

-A 

c.  Thermal  capacitance  -  1,8  x  10  Btu/®F 

Heater  Control:  Solid-state  device  with 

a.  Single-blow  capability  -  100  <  w  <  1000  Ibci/hr, 
aided  ^  -■»  15®F 

CO 

b.  Periodic  capability  -  50  <  w  <  1000  Ibm/hr, 
amplitude  <=«  lO^F 

1,  Large  Motor  -  7.5  Period  <  150  sec 

2.  Small  Motor  -  1.5  <  Period  <  10  sec 

Test  Section:  Plexiglass  drawer  capable  of  accepting  a  test 
core  of  3.25  in.  x  3.25  in.  x  up  to  4.0  in.  long. 

Temperature  Sensors:  Copper-constcuitan  thermocouples 

a.  Inlet  -  1  junction,  28  gage,  ice  reference. 

b.  orifice  -  1  junction,  28  gage,  ice  reference. 

c.  Upstream  -  4  junctions,  0.002  in.  dia.  wire,  ice 
or  mercury  bath  references. 

d.  Downstream  -  4  junctions,  0.002  in.  dia.  wire;  ice 
or  mercury  bath  or  incoming  inlet  air  references. 
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Pressure  Sensors;  Vertical  an'i  inclined  memometers  reading 

inches  of  water.  Taps  are  located  upstream  and  down¬ 
stream  of  the  test  core  and  at  1-D  emd-^-D  from  the 
orifice  plate. 

Flow  Meter;  Sheu:p-edged  ASHE  orifice,  two  plates,  %>.  1  -  1.301 
in.  dia.,  Mo.  2  -  2.201  in.  dia.,  pipe  diameter  -  3.06  in. 

Blower;  1000  Itea/hr'  at  ?0  in.  of  water  static  pressure  xrise. 

Pitot  Tube;  Standeird  pitot-static  tube,  0.1  in.  OD.  A  recession 
in  the  test  rig  wall  is  provided  for  stowage  duriiig  heat 
transfer  tests. 

Circuit  Layout  of  the  Heater  Control 

The  heater  control  consists  of  three  main  sections,  assembled  in 
sequence  and  enclosed  in  an  aluminum  box.  The  sections  are;  front  panel 
Fig.  I-*l,  firing-circuit  boaurd  Fig.  1-2,  amd  power  boaird  Fig,  1-3. 
Schematic  diagrams  for  the  banks  of  capacitors  C2  and  c^  are  ^cwn  in 
Fig.  1-4. 

Operating  instructions  for  the  Heater  Control 

FAST  BESPC»;S£  SIMGI£-BI£IK 

1.  Set  air  flow  rate. 

2.  Adjust  and  C^  according  to  Table  1-2. 

3.  Set  Rj  and  to  0. 

4.  Turn  SH^  to  OM. 

5.  Set  FUNCTION  SWITCH  to  U. 

6.  Adjust  Rj  for  about  20®F  termpature  rise  usii^  airstream 
thermocouples. 

7.  Set  FUNCTION  SWITCH  to  D. 

8.  Adjust  Rg  for  about  5*F  temperature  rise  using  the  then«;'«>uples. 

9.  Switching  FUNCTION  SWITCH  alternately  from  D  to  U  to  D  etc. 
will  produce  the  required  fast  response  St«p-cp  and  Step-Down 
in  airstream  tenperatxjre  entering  the  test  core. 
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NORMAL  £lllGI£-BLOW 


1.  Set  FUNCTION  SWITCH  to  N. 

2.  Switching  SWj^  from  ON  to  OFF  to  ON  etc.  will  produce  a  regl^lar 
response  step-change  in  airstxeam  temperature.  Only  is  used 
to  control  the  upper  tenperature  level. 

PERIODIC 

1.  Set  FUNCTION  St^ITCH  to  P. 

2.  Turn  SWj^  to  OFF. 

3.  Plug  in  the  variable  resistor  with  zero- level  rosistor  at  zero. 

4.  Rotate  cam  to  minimum  resistance  position. 

5.  Turn  to  ON. 

6.  Adjust  zero-level  resistor  for  about  20®f  rise  in  ctirstrean 
temperatu^. 

7.  Rotate  cam  and  note  tanperature  amplitude  fluctuation. 

8.  If  necessary  adjust  the  amplitude  usii^  the  cam  mounted  resistor 
2uid  repeat  steps  4  through  7  until  the  desired  teuperature  wave 
is  obtained. 

9.  The  period  of  oscillation  is  adjusted  using  the  notor  speed 
control. 


w 

Table  1-2 

CAPACITOR  SETTINGS 

^2 

^3 

(Uan/hr) 

Position 

Position 

100 

1.0 

1.5 

130 

1.5 

2.0 

160 

2.0 

2.5 

200 

2.5 

3.0 

250 

3.0 

3.5 

320 

3.5 

4,0 

40C 

4.0 

4.5 

500 

4.5 

5.0 

630 

5.0 

5.0 

800 

5.5 

5.5 

1000 

6,0 

6.0 
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FRCMT  t’ANCL  (BACKSIDE) 


□  WK  TO  flRiND  CIRCUIT  BOARD 
■  WIRE  OMLY  ON  FRONT  (wiL 

Figure  I~1  Backside  of  the  Heater  Control  Front  Panel 

(sh'Twing  position  of  controls  emd  wiring  connections}. 


nRR«-CIRCOIT  BOARD  (BACKSIDE) 


O  TO  FRONT  FAWV 

■  ONLY  ON  FIRKU  CIRCUIT  BOARD 

-  ON  BACKSiOe 

- ON  FRONTSCE 

a  TO  POWER  board 

Figure  1-2  Heater  Control  Firing  Circuit  Board  (shoAfing 

position  of  circuit  elements  and  wiring  connections). 
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APSENDIX  II 


HEATER  THERMAL  CIRCOIT  PARAMETERS  AKD  ECAP  IKPOT 
!Riermal  Resistance 

The  air  film  convective  thermal  resistance  of  the  heater  wires  is 
defined  as 


A  1 
*  (“>„ 

The  total  surface  area  of  the  wires,  of  length  4  heaters  x  56  passes  x 
3.25  inch/pass  and  dla.  0.003  inches,  is 


TT  X  .003  X  4  X  ‘56  X  3. 
:  144 


25  2 

—  =  0.0476  ft 


The  surf ac'i  average  convection  heat  transfer  coefficient  of  the  wire 
h^  can  be  determined  usir.g  data  presented  by  McAdams  [9].  la  the 
reference,  is  plotted  versus  for  air  flowing  normal  to  single 
cylinders.  For  example,  consider  the  air  flow  rate  of  590  IbE^/hr.  The 
Reynolds  ntsober  is  calculated: 


A 

R  \l 


.003  x  590  X  12 


3.25^  X  3600  X  145  x  10~’ 


=  38,5 


Frcaa  the  versus  graph 


3.1 


Air  viscosity  and  thermal  conductivity  used  in  the  and  calculations 
are  taken  from  [iO],  cising  a  film  temperature  of  2G0®F.  Then 


3.1  x  0.0182  X  12 
0.003 


=  226 


Btu 


hr-ft  -»F 
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ais? 


9  =  ' 


^  X  O.C4?£ 


Sta 


Tbe  a2»^  calca.atineis  har^  ssa»  seoes^ia  :cr  xlsar  sasss  cgscai 

tie  carve  iss  rig-  3-2i 

ISsama'!  Ca>anit*y:g 

Valias  used  for  tie  ^svsiczl  ^^i^sztiss  of  alanel  r^atar  aSre 
ass  [131 


=  S<4  Ibs/rtf 


c^  =  0*lj 


32i6-*r 


£he  iotal  casadtaace  da  be 


_  2 

CK5}..  = 
n  4  "  s 


3.14  X  S.C33  X  56  X  3.25  x  544  x  S.i- 


1.77  X  33 


SCap  izspat 

Sbe  3Ca=  f3l  ^zivaied  eircsit:  ior  tie  tbersal  circarb  is  Fig.  ^'- 
is  gitea  in  xi^.,  n-l-  Scdes  ««  pgsibi^  c-Lxrgb  cssectircs  are  assi^ 
SY  the  figsze.  'Sse  SCaP  ir«,:b  fer  tbe  case  ?4*r>6  v  =  5S5  VssJtr, 
considered  in  axsendix  Us,  in  listed  in  ^able  XX-l.  cnrrsEti-sssrr:^ 
'saiaes  on  the  II  c_  rd  are  fcs  the  st^-»  ca»  viin  —  lil,l  pF-  Sbii 
is  discsssed  acre  in  ;<nsedix  XS. 


—  0 


Figure  ii-i  ecat  Representation  of  the  Thermal  Circuit 


Table  ll-l 

ECAP  PROGRAM  INPOT  SAMPIE 


transient  analysis 

H-590  LBN/HR  C?*10.1  HFD 
NCOtl),  C*l.77F-0A 
N(l»0lfR«336. 

C5I  *2«37f2,l0f  1,85,1,64, 1, 46,1. ?2,1 ,71,1,] ^ 

*tii?e*STEP*o 
OUTPUT  INTERVAL* 10 
FINAL  TIMe«o,25 
PRINT,NV,CA 
EXECUTE 


APJEMDIX  III 
STEP-UP  CIRCUIT  ANALYSIS 

Technique  Used 

Ihe  key  elements  In  the  heater  control  circuit  of  Fig.  2.3  can  be 
approximated  by  Pig,  Ili-l; 


22  volts  DC 


Figure  III-I  Approximation  of  the  Heater 
Control  Circuit,  step-Up  Case 

To  analyze  the  circuit  in  Fig.  2.3  a  means  is  needed  to  obtain  the  inter- 
base  voltage  on  the  UJT  as  a  function  of  time.  The  circuit  in  Fig. 
III-l  makes  this  possible  by  simplifying  tlie  calculations.  The  22-volt 
DC  voltage  source  is  equivalent  to  the  voltage  developed  across  zener 
diode  CRg.  interbase  resistance  (i^  =  o)  is  denoted  by  R  .  The  value 
used  for  R^  in  the  design  calculations  was  5910  ohms.  ^ 

The  circuit -equations  written  on  a  nodal  basis  sure 
Node  at  V _ 


■^bb^V^bb  , 

— +  — 5 0 


{III-D 


»nHlHK<H<illl|liMl4i  hUHMi. 


Nofie  at 


V  -V 
BB  2 


+  C, 


“de 


=  0 


Solving  Eqs.  (III-l)  and  (III-2)  for  V  yields 

BB 


(m-2) 


VgB  =20.6-8.7  exp  (-1)0) 


(III-3) 


wl»ere 


The  qucintity  c  is  a  correction  applied  to  the  value  of  to  account 
for  the  chopped  voltage  applied  to  the  firing  circuit,  see  Fig.  2.4.  \ 

chopped  voltage  exists  across  the  firing  circuit  (R^  or  CR^)  because  the 
voltage  collapses  during  conduction  by  the  sCR's.  Fig.  III-2  shows  how 
the  effective  heater  control  output  voltage  varies  with  the  phase  angle 
(triggering  time).  Hie  values  in  the  figtire  are  fxcm  a  sin?>le  heating- 
effect  analysis  of  a  sinusoidal  voltage  wave,  with  a  peak  absolute  value 
of  170  volts.  Figxire  I1I-2  also  defines  the  triggering  time  0^,  used 
in  the  following  definition: 


6  +  Q 

A  s,  initial  Condition  s,  steady-state 

^  2  X  0.00833 


(III-5) 


0  can  be  regarded  as  an  average  portion  of  real  time  during  which 

capacitor  C2  is  being  charged  by  the  blocking  effect  of  the  SCR*s. 

Znterbase  voltage  V  may  now  be  calculated  as  a  function  of  real 
Bo 

time  using  Eqs.  (111-3),  (III-4),  «id  (lll-S).  Trials  are  made  and  the 

results  examined  to  determine  is  a  function  of  air  flow  rate, 

hence  the  variable  notation  in  Fig.  IIl-lv 

Having  obtained  V  (0),  the  firing  voltage  V  may  be  caJ.culated. 

BB  B 

The  relation  used  is  given  by  the  -anufacturer  [7]  for  OJT  2K2646  as 
typically 
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•».  «*e 


III-2 


Effec±ive  Heater  Control  Output 
variation  with  TriggerdLng  Titae 


Vo’ 


igure  III-3 


Heater  Control  Output  Voltage 
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Vp(6)  =  0.65  Vgg(e)  +  0.5 


(III-6) 


Fig.  111^3  providtjs  the  output  voltage  of  the  heater  control  as  a 
function  of  the  firing  voltage  Vp  and  the  steady-otate  resistor  value 
It  is  a  key  one  in  'Qie  circuit  analysis.  Knowing  the  steady-state 
heater  power  required^  the  value  for  ^  found  by  first  using  Fig. 

3.3  to  get  the  steady-state  voltage  outq>ut  required  of  the  heater  control. 
Next,  Fig.  III-3  is  used  to  detesaine  the  value  of  R^  on  the  steady-state 
line  corresponding  to  that  heater  control  voltage* 


SBgf>le  Calculation 

An  intenaediate  air  flow  rate  of  v  -  590  Ibsi/hr  is  chosen  for 
this  exaa^le  with  -  10.1  }iF.  The  heater  power  for  stea<^-state 
operation  is 


*30  = 


=  wc. 


At  = 


590  X  0.24  X 
3600 


20  Btu 

—  =  0.787  ~ 
sec 


OSie  steady-state  voltage  octq>ut  of  the  heater  control  is  then  found  using 
Fig.  3.3: 


H,  Steady-State 


=  64  volts  RMS 


and  the  value  of  R^  froa  Fig.  III-3: 

R^  =  26,000  ohss 

Fie.  m— 3  can  also  be  used  to  find  the  initial  voltage  applied  to  the 
hsateirs.  From  the  figure 


'h.  Initial 


=  112  volts  RMS 


with  a  coirresponding  power  generation  rate  of  2. 37  Btu/sec. 

Big.  III-2  is  used  to  find  the  initial-condition  and  steady-state 
triggering  tises,  9^,  used  in  the  calculation  of  c*  Using  £q.  (Ill— 5) 


0.C0235  +  0.00512 
2  X  0.C0833 


0.45 
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The  value  for  finte  0  is  now  increased  to  0.01  sec.  and  the  series  of 
calculations  beginning  with  V_{e)  is  repeated.  The  results  are  plotted 

DO 

in  Fig.  III-4  for  four  values  of  C^,  covering  the  design  range  of 
interest. 

Each  curve  in  Fig.  III-4  roay  be  used  as  input  for  the  ECAP  program, 
discussed  briefly  in  Section  3.3  and  Appendix  II.  Response  curves, 
calculated  using  ECAP,  are  shown  in  Fig.  III-5  for  the  four  values  of  C^. 
ahe  influence  of  C2  on  the  results  is  very  apparent  and  C2  =  10.1  pF 
wsis  selected  for  use  in  che  circuit. 
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L 


1 


<*l'  S90  Itaa/hr  ,  *  20*F 

ft«»  aos»  «te  .  Bt  •  2S,soon 
K«*  soon 


The  saae  type  of.  analysis  was  done  at  air  flow  rates  of  100  llzs/hr 
and  1000  Ibs/hr.  The  respective  values  for  were  selected  on  the 

basis  of  resulting  closest  a{^roa<di  to  a  step-change  in  airstreaa 
teroerature.  Values  of  to  be  used  at  intermediate  air  floo^  rates 
were  determined  by  interpolation. 


APKNDIX  IV 


STEP-DGSra  CIRCUIT  ANALYSIS 


Figure  XV-l  Ax^ircocusation  of  rte  Heater  Control 
Circuit,  Step-Dcwn  Case 


As  in  tile  step-i^  case,  to  analYze  tj»i  circuit  in  Fig.  2«3  a  aeans  is 

seeoed  to  c^itain  tiie  interbeise  voltage  v  on  tl»  G5T  as  a  fcnction  of 

ss 

tise.  circuit  in  Fig.  IV-1  sisslifies  tne  calculations.  The  22- 

volt  SC  volume  is  equivalent  to  tiie  volt^^e  developed  across  zener 
diode  CR^.  Ynteztase  resistance  (vith  ss  0}  is  denoted  by  The 
value  used  for  is  the  design  calcslatioos  vas  591Q  ohss.  vasiahle 
r^istor  is  used  to  adjust  ^e  stea^^state  i^ltage  fs  S^}  out^^ 
to  the  heaters  is  the  st^i-doifs  aode  of  operati<si. 


. . . . . . . . .  . . . . .  . . . .  ....I  •»..  Mmh,.  . . .  i  .^*>,0  .  .  ji-ft;. 


The  circuit  equations  written  on  a  nodal  basis  are 


Node  at  y 


0-V  V^“V  d(V -V^„) 

BB  .  2  BB  ,  „  2  BB  _  „ 

"rT”  3  ^ 

Q  5 


(IV-X) 


Node  at  V. 


22-^2  ^  W2  ^  ^  ^ 

*~r]"  *  ^3  d8  " 


{IV-2) 


solving  Bqs.  (IV-'l)  and  <IV-2>  for  yields: 
V„(6)  -  [20.6  - \]  o^x-tie) 


■s^  +  •==  +  1 

R-  '  K. 

^  Q  Q  J 


where 


a  is  defined  and  discussed  in  ftppeidix  Ill.  However,  for  the  step-down 
case  the  initial  c<sxdition  corresponds  to  =  0  steady-state 

condition  to  dCT  a-  5®F.  Therefore,  Eq.  (III-5)  becocses 


O  =  0.5  + 


"s,  steady-state 
2  X  0.00B33 


(IV-4) 


As  steady-state  is  approached,  the  exponential  tern  in  Eq.  (IV-3> 
vanishes,  allowing  calculation  of  iC.  Sie  steady-state  heat  generation 


rate  rs 


a  =  wc  r 
*0  P 


ciy-5) 


Fio.  3.3  is  used  to  convert  this  into  heater  voltage.  Fig.  III-3  is 
then  entered  at  the  curve  for  lU  =  40,000  c*ss  to  obtain  the  firii^ 


volt^e  herK:e  V 


B3, steady-state 


using 


V  -0.5 

V  =  — £ - 

BB  0.65 


Fran  Eq.  (IV-3) 


(v 


steady-state 


(IV-6) 


All  the  parameters  in  Sq.  (IV-3)  having  been  evaluated,  the  interfaose 
voltage  V  can  now  be  calculated  as  a  fxmction  of  time.  A  process  very 
similar  to  the  step-up  case  in  Ap^^ndir  III  is  used  to  calctilate  heater- 
wire  temperature.  I5ie  curve  for  R«  =  40,000  dhss  in  Fig.  1II-3  is  used 
for  all  air  flow  rates  in  the  step-down  mode;  the  variable  resistor 
has  been  replaced  by  a  fixed  resistor  Rg. 
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ISHIVATIC^  F<Sl  TES  ".SISiKIgS 

Is  this  a^iendjx,  a  segsssat.  oz  aa  i^alizsd  Isat  exshss^r,  Sis. 

V-1,  is  analyzed  Lo  arrive  art  the  differential  and  intacrsi  easatisss 
vhich  govern  tiie  balk  flaid  tesseratere  aj^  t2s  vail  tanseratere  is  tis 
heat  extdianger.  2fee  c<n:isezvation  of  ei^rsv;  occservatisss  ©f  slsss*  && 
ccndactive  and  convective  2»te  egsaticss,  ai^  the  sgratiac  ef  s^ate  ass 
used  in  tls  asailysis.  fhree  cases  are  cassidered;  laagitadsra^ 

conduction  is  zero  {as  treated  in  [4j),  Sie  veil  lgneit'»->disal  ccri-f^rtim 
is  finite,  and  the  wall  lonsitedinal  ccc^^tlKi  ir  infinite-  in  all 
cases  aie  radial  theraal  candEictivit^  is  considered  as  infinite. 


2.  The  system  is  overall  adiabatic. 

3.  Wall  and  fluid  properties  are  constant. 

4.  The  fluids  are  considered  to  be  lov  velocity  gases  at  constant 
pressure,  and  the  enthalpy  and  internal  thermal  energy  can  be 
treated  as  functions  of  temperature  using  the  specific  heat 
property. 

5.  There  is  no  longitudinal  heat  conduction  in  the  fluid. 

6.  Fluid  thermal  capacitance  is  neglected  relative  to  the  wall 
capacitance  C^.  This  was  considered  in  [4]  and  shown  to  be 
negligible  in  situations  involving  a  gas. 

7.  The  fluid  has  steady-flow. 

Derivation  of  Equations 

The  conservation  of  energy  applied  to  the  solid  wall  segment  of 
unit  depth  normal  to  Fig.  V-1  can  be  stated  by  the  following  equation. 
At  any  instant  of  time  Q 


The  conduction  and  convection  rate  equations  for  g^  and  Aq  respectively 
are: 


w 


g  =  -k  t  _ 
w  9x 


(V-2) 


x+Ax 


Aq  =  J  h(Tj-T^)dx 


{V-3) 


The  energy  storage  term  for  the  wail  is 

x+Ax 

(pet) 


m  =  f 

w  J 


w 


w  99 


dx 


{V-4> 


Eqs.  (V-1)  through  {V-4)  can  be  combined  and  cifter  applying  the  theorem 
of  the  mean  and  passing  to  the  limit  (Ax  -»  0),  the  following  equation  is 
obtained; 
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(V-5) 


i  .  aT  V  a-r, 

-’'I:  -w  “  ° 

The  conservation  of  energy  applied  to  the  fluid  system  in  Fig.  V-1 
can  be  stated  by  the  following  equation.  At  any  instant  of  time  6 

(wi)x^^  -  (wi)x  +  Aq  =  0  (V-6) 

Eqs.  (v-3)  and  (V-6)  can  be  combined  and  after  applying  the  theorem  of 
the  mean  eind  passing  to  a  limit  (Ax-^  0),  the  following  equation  is 
obtained: 

-^(wi)  +  h(T  -T  )  =  0  (V-7) 

O^v  X  V? 


For  a  perfect  gas 


c  ^ 
ax  p  ax 


(V-8) 


since  w  is  a  constant,  Eqs.  (V-7)  and  (V-8)  are  combined  and  become 


aTf 

wc  -r —  +  h(T^-T  )  =  0 
p  ax  f  w 


(V-9) 


The  parameters  C^,  C^,  ^d  R  are  now  introduced  and  are  defined  by 


C  -  ^  wc 
f  =  p 


rA^ 
^  =  hA 


Eqs.  (V-5)  and  (V-9)  become 

it)  ~ 


Tf  -  T  .  bT 
^  ”  +  c  _ -  =  0 

K  w  ae 


(V-IO) 


aTf 

ax  RL 


(V-11) 
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I4ie  above  Eqs.  (V-IO)  and  (V-11)  apply  even  for  variable  A^,  R,  C^,  and 
C^.  These  will  be  considered  constemt  and  the  following  definit:ions 
are  introduced: 


*  AX 
=  L 


0* 

RC, 


w 


RC  ^  ^ 

=  N. 
tu 


k  A 

i  A  _«ji 
^  =  rf; 


{V-12) 


m  _r,'« 

rn*  A  f  ‘ro 

f  =  AP 


m  _m 

*  A  w  m 

T  =  . 

w  ZJtr 


2  =  * 


where  T  is  the  mean  temperature  about  which  the  system  oscillates 
and  At  is  the  amplitude  of  T,  at  the  inlet.  T,  and  T  are  the  non« 
dimensional  temperatures  and  9  is  the  nondimensional  time.  Bie  non- 
dimensional  space  coordinate  is  2.  Then  Eqs.  (V-IO)  and  (V-11)  become: 


*  2  * 

.  *  3T  d  T 

if  *  W  W 

®f  '■^w—'^^u  — 
36  3Z 


(V-13) 


*  * 

T-  -  T  =  -  -—r 
f  w  92 


(V-14) 


Solution  For  Wall  Longitudinal  Conduction  Zero  (X  =  0) 

/  h  A  V 
/a  W  W  I 

when  the  wall  longitudinal  conduction  parameter  X^=  ) 

approaches  zero,  the  last  term  in  Eq.  (V-13)  becomes  negligible  and  the 
equations  of  interest  are 


*  *  3T 

it  *  W 

«x»  -  X  =  — 

^  36 


{V-15) 


^ 

’’f  "  "^w  ■'  92 


(V-14) 
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The  boundary  condition  is 


T ,  ( 0, 9  )  =  sin2n  =  sin  ~ 


{V-16) 


with  the  definition 


8. 


•  ^  ^  ~o 
^  =  2tt  Rc„ 


{V-17) 


The  parameter  ^  is  a  nondimensiona 1  period  of  oscillation  and  RC^  is 


the  heat  exchanger  wall  time  constant*  The  solutions  of  interest  are 
the  particular  solutions  which  result  after  all  the  initial  condition 
effects  have  damped  out.  The  solutions  cis  given  in  [4]  au:e 


(V-18) 


T*  =  »  ■  exp/ — I  sin  - ^ — -  Z  -  a]l 

«  JTTp  u  ,  +  If 


(V-19) 


where  B  ^  tan”^  -i 
w 


(V-20) 


For  the  fluid: 


Amplitude  Attenuation  =  =  exp^ - ^72 ) 


{V-21) 


and 


thase  Shift  = 


Yf 


(V-22) 


The  amplitude  octer.uation  eind  the  phase  shift  of  the  fluid  at  x  =  1 


or  Z  =  are  plotted  in  Figs.  V-2  through  V-9. 


For  the  wall: 


9  /  *2*  \  T  v>* 

Aciplitude  Attenuation  =  d  -  /■  exp  {  ■  ■■  i  =*  >  ■  h 

“  '1+4^'  CTr  * 


>«  . . . 


Phase  Shift  for  the  Fluid  Figure  V“3  Phase  Shift  for  the  Fluid 

Temperature  Response,  Temperature  Response, 

X*  “  1  A  »  1. 0.  X*  »  1,  ’If  ®  2,0, 
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Phase  Shift  for  the  Fluid  Figure  V-5  Phas  >  Shift  for  the  Fluid 

Temperature  Response,  Temperature  Response, 


CV-24) 


iaDi.e  v*-i 


CCNSTANT  cilN  TH6  SOLUTION  FOR  TH6  PERIODIC  TFCHNiQUE 
NIOj  finite  mall  LuNGlTtK)lN4L  CunOUCTIuN 


K3t=:rs 

n — 

H 

II 

H  O 

H 

II 

0.001 

0.005 

O.OlO 

0.050 

O.iOO 

0.500 

l.OOO 

0,.2 

0.9t>l5 

0.9616 

0.9617 

0.9618 

0.9629 

0.  9642 

0. 9767 

0.  9974 

C.4 

0.8621 

0.8621 

0.6621 

0.8621 

0.8623 

0.8625 

0.8533 

0.7989 

0.6 

0.7353 

0,7352 

0.7347 

0. 7342 

0.  7296 

0. 7232 

0.6506 

0,5498 

0.8 

0.6098 

0.6095 

0.6087 

0.6076 

0.5988 

0.5875 

0-5000 

0.4242 

l.O 

C.5000 

0.4998 

0.4983 

0.4975 

0.4877 

0.4758 

0.4003 

0.3457 

1.2 

G.4096 

0.4096 

0.4087 

0.4076 

0.3989 

9. 3889 

0.3313 

0.291? 

1.4 

0.3378 

0.3377 

0.3370 

0.3361 

0.3294 

0.321B 

0.2799 

0.2506 

1.6 

0.2809 

0.2808 

C.2803 

0.2797 

0.2750 

0.2697 

0.2404 

0,2191 

1.8 

0.2358 

0.2358 

0.2355 

0.2351 

0.2321 

0.2287 

0.2091 

0.1939 

2.0 

C.2000 

0.2000 

0.  ! 998 

0.1996 

0. 1 979 

0.1959 

0.1833 

0.1733 

2,2 

0.1712 

0.1712 

0.1711 

0.1710 

0. 1 703 

0.  lo94 

0.1623 

0.1560 

2.4 

0. 1479 

0.1479 

0.1479 

0,1479 

0.1478 

0.1477 

0.1453 

D.1413 

2.6 

0.1289 

0.1289 

0. 1289 

0. 1290 

0.1294 

0. 1298 

0.1305 

0.1288 

2.8 

0.1l3i 

0.1131 

0.1132 

0.1133 

0.1140 

0.1148 

0.1178 

9. 1179 

3.G 

O.lOQO 

0.1000 

O.lOOl 

0.1002 

0.1012 

0.1022 

0,1069 

0.1083 

3.2 

0.0890 

0.0800 

0.0891 

0.0892 

0.0903 

0.0915 

0. 0974 

0. 0999 

3.4 

0.0796 

0.0796 

0.0798 

0.0799 

0.0811 

0.0824 

D.0891 

0.0925 

3.6 

0.0716 

0.0717 

0.0718 

0.0719 

0.0731 

0.0745 

0. 0818 

0-0858 

3.8 

0.0648 

0.0648 

0.0649 

0.0651 

0.0663 

0.0676 

0-0753 

0. 0799 

4.0 

0.0588 

C.05S9 

0.0590 

0.  0591 

0.0603 

0.0617 

3.0695 

0.0745 

4.2 

0 .0536 

C.0537 

0.0538 

0.0539 

0.0551 

0.  0565 

0.0644 

0.0697 

4.4 

0.0491 

0.0491 

0.0493 

0.0494 

0.0505 

0.0519 

0.0598 

0.0653 

4.6 

0.0451 

0.0452 

0.0453 

0. 0454 

0.0465 

0.  0478 

0.0556 

0.0oi3 

4.8 

0.0416 

0.0416 

0.0417 

0.0419 

0.0429 

0.0442 

0.0519 

0. 0576 

S.O 

0.0385 

0.0385 

0.0386 

0.0387 

0.0397 

0. 0409 

0.0435 

0.0543 

5.2 

0.0357 

0.0357 

0.0358 

0.0359 

0.0369 

0,0330 

0.0454 

0. 0512 

5.4 

0.0332 

0.0332 

0. 0333 

0.0334 

0.0343 

0.0354 

0 .0426 

0.0484 

5.6 

0.0309 

0.0309 

0.0310 

0.0311 

0.0320 

0.  0331 

0.0401 

0,0458 

5.8 

0.0289 

0.0289 

0.0290 

0.0291 

0.0299 

0.0319 

0.0377 

0.  04  34 

6.C 

0.0270 

0.0270 

0.0271 

0.0272 

0.0280 

0. 0290 

0,0356 

0.0412 

6.2 

0.0254 

0.0254 

0.0255 

0.0256 

0.02O3 

0.0273 

0,0336 

0.  0391 

6.4 

0.0238 

0.0239 

0.0239 

0.0240 

0.0248 

0.0257 

0.0318 

0.0372 

6,6 

0,0224 

0.0225 

0.0225 

0.0226 

0.0233 

0. 0242 

0.0301 

0.0354 

6.8 

0.0212 

0.02 12 

0.0213 

0.0213 

0.0220 

0.922R 

0.0235 

0.0337 

7.0 

0.0200 

0.0200 

0.0201 

9.0202 

0.  0208 

0. 0216 

0.0271 

0.0322 

7.2 

0.0189 

0.0189 

0i0l90 

0.0191 

0.0197 

0.0205 

0.0258 

0,0307 

7.4 

0.0179 

0.0179 

0.01 SO 

0.0181 

0.0187 

0.0194 

0.0245 

0,0294 

7.6 

0.0170 

0.0170 

0.0171 

0.0172 

0.0177 

0.0184 

0.0234 

0.  0281 

7.8 

0.0162 

0.0162 

0.0162 

0.0163 

0,0169 

0.0175 

0.0223 

0,0269 

8.0 

0.0154 

0.0154 

0.0155 

0.0155 

O.OISO 

0-0167 

0.0213 

0.0258 

3.2 

0.0147 

0.0147 

0.0147 

0.0148 

0.0153 

0.0159 

0.0203 

0.0247 

8.4 

0.0140 

0.0140 

0w0140 

0.0141 

0.0146 

0.0152 

0.0195 

0.0737 

8.6 

0,0133 

0.0134 

0.0134 

0.0135 

0.0139 

0.0145 

0. 0136 

0.0228 

8.8 

0.0127 

0.0128 

0.012S 

0.0129 

0.0133 

0-0139 

0.0179 

0.0219 

9.0 

0.0122 

0.0122 

0.0122 

0. 0123 

0.0127 

0. 0133 

0.0171 

0.0211 

9.2 

0.0117 

0.0117 

0.0117 

O.OllB 

0.0122 

0.0127 

0.0164 

0,0203 

9.4 

0.0112 

0.0112 

0.0112 

0.0113 

0.0117 

0.0122 

3.0153 

0.0195 

9.6 

0.0 10 y 

^.0107 

0.0108 

0.0108 

0,0112 

0.0117 

0,0152 

0.  01  S3 

9,8 

0.0103 

0.0103 

0.0104 

0.0104 

0,0183 

0.0112 

9.0146 

0.0132 

n 

II 

II 

II 

II 

II 

II 

II 

M 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

II 

u. 

85 


86 


Newton-Raphson  eethod  [16]  to  fins?  the  roots  and  ^ 

used  to  calculate  and  for  intermediate  valoes  of  f  asd  is^  . 

not  provided  by  the  tables.  It  is  entered  aeans  of  the  fortran 
statement,  see  JU^Mndix  VII.  For  the  fluid 


and 


Khen  1  =  0 


^^litude  Attenuation  =  D- 

X 


e»{-c^S} 


Ifcase  Siift  =  =  c^Z 


and 


1  +  i 


1  ■*'  t" 


CV-2S} 


Solution  For  sail  loagitcdic*!  Ccr.doctiaa  Infinite 

As  the  vail  Icsmgitud  Inal  contfoctica  paraseter  X  becones  very- 
large  (X -*■«»),  the  te-j£perat-2re  T  of  the.  wall  becxssss  only  a  fsKrtim 

m 

of  0,  The  te^^tsre  of'  the  fXiiJ,  ressaim.  a  fur^rtire  cf  both 
tttre  9  and  distance  x.  Then  tv-S)  and  fV-9)  br.-oEC 

d? 

+  fpct)^  ^  *  Q  {¥-31J 


VC 


P 


+  h  CT^-T  }  =  0 

■  X  V 


nr-22i 


•She  canservation  of  energy  allied  to  the  ^lid  vaXI  ssgassz  cf  =nit 
depth  nocsal  to  Fig.  v-1,  can  be  stated  by  the  foilcadng  eqsatice: 


I. 

-h  f  Idx  +  Lisct) 

J  *.  w 
o 


V  ^ 


331 


The  conservation  of  energy  applied  to  the  flsic  systea  iii  Fig.  ^1  css 
be  stated: 


dx  -T-  h 


/Cr,-T  jdx  =  G 

i.  V 

o 


8? 


. . . 


introducjrig  some  of  the  nomenclature  of  Eqs.  (V-12),  Eqs.  (V-23)  and 


(V-34)  become; 


“tu  * 

1  r  *  *  <3^v, 

■  7r~  /  T,  dZ  +  T  +  -—X  =  0 

j  £  »  le* 

o 


(V-35) 


f*""  afTf  1 


•  dZ  -  T  =  0  (V-36) 

I  w 


The  solutions  arr. 


♦  V  e  z  * 

T*.exp(-e3z)sil.(^)+^cos(^)  (V-37) 


where 


®3  N 


(V-33) 


(V-39) 


(V-40) 


* 

The  phase  shift  of  tiie  fluid  for  \  =  co  at  x  =1  is  presented  in 
Figs.  V-2  through  V~5.  The  amplitude  attenuation  of  the  fluid  for 

■k 

\  ~  a>  at  X  =1  is  presented  in  Figs,  V-6  through  V-9. 

Figs.  V-10,  V-11,  and  V-12  show  what  is  happening  to  both  the  fluid 

•  *  S* 

ten^rature  T_,  and  the  wall  temperature  T  ,  as  functions  of  time 

space  X  ,  and  wall  longitudinal  conduction  Fig.  v-10  shows  the  time- 

* 

varying  temperatures  at  the  heat  exchanger  matrix  exit  (x  =1).  The 


38 


r 

0.09 

QOS 

ao7 

b 

0.06 

r 

0.05 

Q04 

003 

i 

QQ2 

f 

00) 

1 

0 

i 

-00) 

-002 

-003 

h.  ^  ^ 

dyfj/  c»  zrr  0/6t  > 


Figure  V-12  Sample  of  the  Time  Veuryiug  Ten^rature  Difference 
(Tf  -  tJ)  at  the  test  matrix  exit,  x*  =  1.  (T*--? 
plot  is  phase-shifted  in  time  to  appecu:  as  a 
sinusoidal  curve  beginning  at  zero.) 


changes  in  ai[j)litude  attenuation  and  phase  shift  due  to  longitudinal 

conduction  effects  can  be  seen*  Fig.  V-11  shows  the  wall  and  fluid 

* 

temperatures  along  the  longitudinal  eucis  of  the  core  at  time  0=0. 
Again  the  effects  of  longitudinal  conduction  can  be  seen,  however  in  this 
case  the  amplitude  attenuation  change  and  the  phase  shift  chemge  cannot 
be  read  directly.  Pig.  V-12  shows  the  temperature  difference  between 
the  fluid  and  the  wall  for  the  (X  =  0)  case.  The  difference  is  a  direct 
measure  of  heat  flux  to  2uid  from  the  wall.  Because  it  is  small  and 
creates  an  illusion,  it  is  difficult  to  read  this  difference  from  the 
X  =  0  curve  in  Fig.  V-10.  A  dotted  line  representing  the  constant  heat 
flux  boundeury  condition  is  superimposed  on  the  figure.  It  shows  that 
the  periodic  technique  tends  toward  the  constant  heat  flux  boundary 
condition.  Results  of  testing  are  therefore  not  expected  to  agree  fully 
with  theoretical  predictions  using  the  constant  heat  flux  boundary  con¬ 
dition. 
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APPENDIX  VI 


ERROR  ANALYSIS  FOR  TOE  PERIODIC  TECHNIQUE 


The  periodic  treinsieiit  technique  requires  accurate  measurement  of 

temperature  variation  upstream  and  downstream  of  the  core.  Matrix  j^ysi- 

cal  properties  such  as  density  p  ,  specific  heat  c  ,  cuad  thermal 

w  w 

conductivity  must  also  be  known.  The  t^jective  of  this  appendix 
is  to  show  how  uncertainties  in  the  temperature  measurements  and  in  the 
core  physical  properties  may  aiffect  the  heat  transfer  test  results  (the 
j — N  curves) . 

X\ 

Fluid  Temperature  Measurement,  Sinusoidal  Case 

The  solution  for  the  fluid  temperature  when  wall  longitudinal  con- 

* 

duction  is  zefo  (X  =  0)  is  given  by  Eq.  {V-18) .  At  x  =  1  this  equa¬ 
tion  becomes 


Amplitude  Attenuation  Method  -  The  amplitude  attenuation  method  compares 
the  amplitude  of  the  temperature  wave  downstream  of  the  test  core  with 

it 

the  input  wave  upstream  of  the  core.  At  x  =  1  the  amplitude  attenua¬ 
tion  from  Eq.  (VI-1)  is 

-N 

D,  =  exp  ( - 1-)  (IV-2) 

'1  +  ^  ' 


Experimentally  this  is  arrived  at  by  recording  tlie  upstream  and  down¬ 
stream  temperature  waves  and  selecting  the  maximum  and  minimum  values  on 
the  waves.  If  the  outlet  amplitude,  called  ^^2'  error,  the 

following  result  is  obtained  from  Eq.  (VI-2) : 


AT 


f2 


At 


{VI-3) 
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and 


dh/h 

f2  o 


2  ^ 
It  -w 


(VI-4) 


Note  the  singularity  in  Eg.  {VI-4),  for  ^  =  1  the  amplification  factor 
goes  to  eo. 


Phase  Shift  Method  -  The  phase  shift  method  compares  the  phase  of  the 

tenperature  wave  downstream  of  the  test  core  with  the  input  wave  upstream 

* 

of  the  core.  At  x  =1,  the  phase  shift  from  Eg.  (VI-1)  is 


■  1 1 


N^ 


(VI-5) 


E}q>eriinentally  the  phase  shift  is  arrived  at  by  determining  what  portion 
of  the  input  wave  length  (considered  as  2tT  radians)  has  been  shifted  by 
passage  of  the  fluid  through  the  test  core.  If  the  outlet  temperature, 
called  Tj2'  hhe  following  result  is  obtained  from  2g. 

(VI- 1) ; 


/vr 

o 


(VI-6) 


Noting  that  the  phase  shift  is  measured  at 


1  +  t 


■  tu 


(VI-7) 


it  follows  that 


dh/h  ^  JL_  r_yL+lL.1 

dTfz/^-o'Df  |-2N^^^J 


(VI-8) 


For  a  given  N^^,  one  may  evaluate  the  right-hand  sides  of  Egs. 

(VI-4)  and  (VI-8)  for  various  values  of  ?.  The  best  experimental  value 
of  f  (with  least  resulting  xmcertainty)  and  the  best  method  (.amplitude 
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attenuation  ox  phase  shift)  stay  be  selected  and  used  in  testing.  Fig. 
7.3  is  an  operationally  convenient  plot  showing  the  values  of  t 
tlie  ajethod  to  use  for  values  of  0.2  <  <  50.  This  figure  is  a  plot 

of  i  versus  and  the  be.'.  t  test  method  that  will  yield  the  minimum 

uncertainty  for  the  heat  trans.cer  coefficient  is  noted  on  the  graph. 
The  magnitude  of  this  uncertainty  is  about  3  percent  on  h  for  a 
dr  magnitude  of  1  percent. 

O 

Test  Core  ihysical  Properties 

The  periodic  transient  technique  requires  that  the  siatrix  £^ysical 

pecific  heat  c^,  and  thermal  conduc- 
Tha  technique  calculates  the  parameter  f  from 
experimentally  obtained  and  using  Eqs.  {VIl-17)  and  {VII- 18) . 


properties  such  as  density  p^, 
tivity  Jc^  be  known 


Density  of  the  Core  Solid  Material  -  The  laboratory  technique  used  for 
calculating  the  core  porosity  p  is  to  measure  the  core  mass  M,  the 
frontal  etreu 
<VI-9) : 

P  *  T- 


the  flow  leitgth  L.  and  use  the  following  Eq, 


*  1  - 


M 


fr 


’fr 


^9. 


(VI-9) 


The  sensitivity  of  the  N  .  to  uncertainty  in  p  may  be  found  using 

O  V  ^ 


where 


«  h  _  c  * 

**St  =  Gc  A  ^  “ 
P 


C  p 
A  “  Ob 


(Vl-iO) 


(VI-11) 


For  the  revised  Bay  ley  Model  [15] 


A  r 
c  _  vp 

A  “ 


(VI-12) 
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where  1]  is  a  constant  based  on  observed  or  neasured  cell  characteris¬ 
tics  such  as  cell  count  or  cell  ^ape.  Cceablning  Egs.  {VZ-9),  {VI- 10), 
and  {VI-12)  »nd  differentiating  yields 

^t  _  ^  _  1  -  P  (IV-13) 

«st 


The  sensitivity  of  to  uncertadnti“E 


4r.  G 
n 

M 


4LG 

P 


in  nay  be  found  using 

A 

{IV-14) 


Combining  Eqs.  {VI-12)  euid  {VI-14)  and  differentiatii^  yields 


(IV-15) 


For  the  ceise  of  fully  de^«loped  laminar  flow  in  tubes,  with  the 
constemt  heat  flux  or  constant  wall  tec^rature  boundary  condition,  and 
a  fixed  geoBetry,  the  Kusselt  number  is  constant  and 


‘Shi 


St  R  Pr 


(IV-IS) 


The  relationship  between  j 


and 


has  the  form 


j  = 


(VI-17) 


where  E  depends  on  the  cross  section  geometry  and  the  boundary 
condition,  and 


dN 

R 


(VI-i8) 


C«;^)iniiw  Eqs.  (VI-13),  (VI-IS),  and  (VI- 18)  produces 


E 


0 


(Vl-19) 


Ifiis  Beans  the  characteristic  is  pot  sensitive  to  uncertainties 

in  even  though  j  emd  if  considered  individually  are. 

Specific  Heat  of  the  Core  Solid  Material  -  The  core  specific  heat  c^ 
will  affect  only  the  heat  transfer  results,  j,  of  the  periodic  tech¬ 
nique.  Elow  friction  results  and  Reynolds  nacher  are  not  affected. 
Differentiatii^  Eq.  (VX-IO)  yields 


dH  * 

.St  ^  ^ 

St  ^  SD 


Fr<xa  the  definition  of 


*  c 
OJ  w 


(VI-20> 


(VI-2i) 


For  the  anq)litude  attenuation  method,  the  equation  used  for 
ccstputing  f  is  Eq.  (VII- 1?).  Csing  Eq.  (VII- 17)  in  Bq.  (VI- 20)  yields: 


■  ■  24^ 


(Ihe  plus  sign  is  used  for  t  >  1,  tiw  ninus  sign  for  t  <  i*) 

For  the  {hase  shj.ft  laethod,  if  is  c.oi:5>uted  using  Eq.  (Vll-18). 
Combining  Eq.  <VII-18)  and  Eq.  (VI-20)  produced: 


1  /x; 

w  ^  V 1  +  t' 


{VI-23) 
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The  sensitivities  in  Eqs.  {VI-22)  and  (vi-23)  are  plotted  in  Fig.  vi-l 
for  0  <  f  <  10,  X  =  0,  and  no  temperature  measurement  error. 

For  the  meucimun  slope  method  [12]  with  X  =  0 


^St'^St 

dc  /c 
w  w 


2.2 


for  the  range  >  10. 

ISiermal  Conductivity  of  the  Core  Solid  Haterial  -  solid  wall  thermal 
condiictivity  affects  only  the  heat  tramsfer  results,  through  the  wall 
longitudinal  conduction  parameter  X*  For  the  periodic  technique,  this 
affect  must  be  taken  into  consideration  if  0.01  {For 

XK^^=  0.01,  the  influence  on  amplitude  attenuation  auid  pliase  shift  is 
less  than  1  percent) .  Generally  the  magnitude  of  the  wall  longitudinal 
conduction  effect  on  heat  transfer  results  is  less  than  10  percent. 
Therefore  an  uncertainty  in  solid  thermal  conductivity  k^  of  10  percent, 
hence  in  X  also,  would  produce  only  about  a  1  percent  uncertainty  in 
heat  transfer  results. 

Consider  for  example  the  actual  test  point  for  Core  510  which  had 
the  largest  XN^^=  0.091.  If  the  effect  of  longitudinal  conduction  was 
canpleteiy  ignored,  the  heat-  transfer  result  would  be  6. 9  percent  in 
error,  under  the  same  test  conditions,  liie  maximum  slope  method  would 
have  been  7.1  percent  in  error  if  the  longitudinal  conduction  effect  was 
cQB^letely  ignored.  Therefore  the  periodic  tichnique  is  about  as  sensi¬ 
tive  as  the  maximum  slope  single-blow  techniqu-  *"0  wall  longitadinal 
conduction  effects. 
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A?SE3S)IX  VII 


DAia  asDocnos 


Single-Blow  iteehnlque 

iSie  coaputer  prograa  used  is  as  de^rribed  by  iSieeler  in  [  5j.  Swsi- 
cal  pr^erties,  suc'  as  specific  tiiexsal  coadisctivityj  jstm  fensity 

of  the  core  solid  saterial  were  t:aken  frees  [l5j  are  ffer  mrtn, 
teisperature  conditions) 


=  141  Ib^ft^ 
c^  =  0.200  3t5j/{lhs  *F) 
=  0.42  Btu/ {hr  St  *!■) 

m 


The  physical  properties  of  air  used  in  the  data  redcstica  srucraa  are 
froc  Bureau  of  Standards  Sources  TlOj  arv^  the  hsiidity  ccn^cticss  ar** 
nade  in  accordance  with  the  scthod  described  by  Shah  Ci2j. 

^e  naxLcas  slope  nethod  is  ised  in  this  progracs.  in  r-^v-iniTtf^’-sdry. 
fon?.  the  naxisaa  slope  is 


Kax.  Slope  =  KS  = 


C  X  36C0 

V _ 

X  {ST3U3} 


DTAO  is  tha  tine  (sec.)  required  for  the  nc-ndlnenSicnal  teacerature  to 
change  its  full  value,  1.0,  at  the  naxisass  rets  of  ^snge.  The  erersge 
of  four  runs  at  a  single  flew  rats  ;’two  step-=^  and  two  step-dcsc'l  is 
used.  The  longitudinal  conduction  perssetsr  \  is  also  required  to 
establi^  K  .  A  graphical  and  tabular  presentation  of  H  as  e  fcrcti^ 

v»V» 

of  Max.  Slope  and  X  is  given  by  fl2].  The  ccs«tar  piagract  iocorporatss 
this  in  the  fora  of  a  "look-up"  table. 

jScrssenclature  for  the  coaputed  data  is  presented  in  Table  vXI-1. 
Single-blow  heat  transfer  test  result  listings  are  provided  for  cores 
SOSA,  5053,  505C,  and  510  in  Table  VIl-2  for  both  a:^  the 

un£ii;^d  cases,  “nje  flow  friction  results  are  virtuallv  iisnticai  nor 


the  aided  and  unaided  cases.  Flow  friction  results  for  core  505A  may 
be  found  in  [15].  Flow  friction  results  for  core  510  are  included  with 
the  periodic  resulus.  in  this  appendix. 


Table  VII- 1 

NOMENCIATURE  FOR  SINGLE-BICW  HEAT  TRANSFER  RESULTS 


W 

G 

TEA 

DTi'U 


LAMBDA 


MAX  SLOPE 

NTU 

NSl 

J 

NR 


Flowrate,  Ibm/hr 

2 

Mass  Velocity,  lbm/(hr  ft  ) 

Average  core  inlet  temperature,  "F 

®ie  time  required  for  the  nondimensional 
temperature  to  change  its  full  value,  1.0, 
at  the  maximum  rate  of  change,  sec 

X,  heat  transfer  longitudinal  conduction 
pareimeter,  dimensionless 

Meucimum  slope,  dimensionless 

N.  ,  dimensionless 
XXi 

Ng^,  Stanton  Number,  dimension? -as  s 

2/3 

Colburn  j  factor  ,  dimensionless 

5  w  Pi 

Reynolds  Number 
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Table  VIl-2  cc.it'd, 
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DATA  REDUCTION  -  cont’d 


Periodic  Technique  -  Heat  Transfer  Calculations 

The  measured  inlet  and  outlet  temperature  variations  (x*  =  o  and  1) 
of  the  fluid  are  processed  using  Eqs*  (V-29)  and  (V-30) 


Df  = 


{V-29) 


Yf  - 

to  determine  the  heat  transfer  performance.  The  terms  and  are 

listed  in  Tables  v-1  a  1  v-2  as  functions  of  and  • 

Since  the  inlet  and  outlet  temperature  waves  are  not  in  general  pure 
sine  waves,  the  analysis  of  Appendix  v  is  modified  using  Fourier  Series 
to  account  for  this.  If  only  the  first  harmonic  is  considered,  a  general 

periodic  temperature  variation  at  the  test  core  inlet  can  be  represented 
by; 

®o  2Tr  5., 

"^(O/Q)  =  -^  +  A  sing—  0  +  E  vos~  0  (VII-1) 

”o  ®o 

where 

®o 

A  =  ^  f  T  (0,0)  sin|li0  de  (Vii-2) 

o  J  Oq 


2  r 


(0,0)  COS'^0  dg 
o 


1  / 

r"  T  (0,0)0  0 

J  ^ 


{VII-3) 


{VII-4) 
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The  measured  periodic  response  of  the  fluid  (x*  =  1)  can  be 
in  like  manner  by* 


where 


'*  o  2tt  2Tr 

T^'x  =  1,  8)  =  -^  +  a  siri— 0  +  b  cos— 0 

®o  °o 


n;ii-5) 


T,{x  =1,9)  sin^8d0 
®o 


(VII-6) 


b  = 


r  * 

/  Tf(x  = 
'O 


=  1,0)  cos— Od0 
®o 


(VlI-7) 


m  2  0^  J 


T^(x  =1,0)  d0 


(VII-8) 


Eqs.  (VII-4)  and  (VII-8)  serve  as  an  e:q)erimental  check  on  the  adiabatic 

*  ★ 

behavior  of  the  test  section  from  x  =  0  to  x  =  1. 

iSie  above  analysis  is  similar  to  that  in  [4]  where  tlie  higher 
harmonics  were  also  considered.  However,  the  higher  harmonics  have  no 
value  in  the  heat  transfer  calculations,  so  they  are  deleted  from  this 
analysis. 

The  Fourier  Coefficients  in  Eqs.  {VII-2),  (VJI  ,  {VII-6),  and 

(ViI-7)  can  be  numerically  determined  using  the  Trapezoidal  Rule.  The 

* 

temperatures,  T^(0,9)  and  T^(x  =  1,6),  roust  be  known  simultaneously 
as  functions  of  time  for  at  least  one  period  9  f  this  rriay  be  done 
graphically-  Each  period  9^  of  these  two  wave  functions  is  divided 
into  (K-1)  equal  intervals  with  K  data  points-  The  first  data  point 
is  considered  as  at  9/9^  =  0  and  the  last  as  at  9/9^^  =  !•  However 
the  first  point  may  occur  at  any  point  in  real  time  on  the  wave  function 
records,  keeping  it  simultaneous  for  both  the  upstream  and  downstream 
records-  The  two  temperature  wave  records  are  used  in  Eqs.  (VII-2), 
{vil-3),  (VII-6),  and  (VII-7)  to  obtain  the  Fourier  Coefficients.  This 
is  done  numerically  using  the  following  approximation  (Trapezoidal  Rule) 
with  B  as  cin  example; 
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where 


T^(0,9}  «  T^(o,^(i-l))  =  0)  i  -  1,2, 

»  1  ^kl 

Vi  .  Tji  cos[2tl(^)]  i.  i,.- . K 


K 

(VlI-9) 

(VII-10) 

(\rti~ii) 


other  Fourier  Coefficients  (A,  a,  and  b)  can  be  calculated  in  like 
manner  by  using  the  appropriate  temperature  record  in  Eq.  (Vll-9)  and 
the  associated  trig  function  in  Eg.  (VII-11) . 

The  number  of  data  points,  K,  typically  ranged  from  30  to  50  for 
the  results  included  in  this  appendix  on  Core  510.  This  makes  the 
Trapezoidal  Rule  reasonably  accurate.  The  number  K  can  be  either  even 
or  odd. 

iSie  temperature  records  were  obtained  with  one  ixonsywell  Electronik 
19  Single-Channel  Lab  Recorder,  To  obtain  the  upstream  and  downstream 
ten5)erature  records,  the  recorder  wais  connected  to  the  upstream  thermo¬ 
couples  for  at  least  one  period  0^  and  then  quickly  switched  to  the 
downstream  set  for  at  least  one  period.  Using  the  sv/itching  »-s;  the 

steurting  point  for  the  downstream  record  and  then  "backtracking"  one 
period  for  the  upstiream  steirting  point,  two  simultaneous  wave  traces 
were  obtained. 

The  rouriar  Coefficients  are  used  to  determine  the  experimental 
amplitude  attenuation  and  phase  shift  as  follows  [4]: 

r”o  2 

Amplitude  Attenuation  =  D_  ^  =-  x'VTI-12) 

Shift  =  (VII-13> 

b3 
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These  relationships  are  combined  with  Eqs.  (V-29)  and  (v-30)  to  obtain 


2  .2 
a  +  b 


177 


(VII-14) 


'1  t. 

-1  /  b  “  5 


bB 


(VII«15) 


3i-r.ce  e,  and  are  complicated  functions  of  i  and  XN  and 
12  ^  tu 


«  2n  A  .  * 


(VII- 16) 


the  above  eq'iations  ceinnot  be  solved  explicitly.  For  the  case  where 

X  =  0,  using  the  amplitude  attenuation  obtained  experimentally  and 
2 

=  1/(1  +  ^  ),  Eq.  (V2I-14)  can  be  solved  for  i; 


*  ( 

r  * 

^  ) 

CD  /2  ^  1 

»  /2 

-  1 

J  ^  1 

^  ( 

\D.  /  1 

'  I,  expr' 

L^iD  / 

'“f,  expr  '* 

1/2 


(VII-17) 


(The  plus  sign  is  used  for  f  >  1,  the  minus  sign  for  t  <  !•)  Using 

2 

the  E^iase  shift  obtained  experimentally  and  62  =  y/d  +  ?  )»  Uq. 
(Vti-15)  can  be  solved  for  y; 


t  =  -F 


- 

^z,exj?r  J 


J./2 


{VII-i8> 


A  heat  transfer  cerformance  narise&tec,  N  , 

~  tu 


can  now  be  calculated  using 


Thesa  values  for  J  and  are  only  fully  applicable  if  m  „  c 
However,  they  ere  the  initial  values  to  use  when  entering  ,Sles  V-1  ».a 
V-2  to  obtain  approrlnate  values  for  Sj  and  e^.  o.e  experi™,ntal 
enplltude  attenuation  and  phase  shift  y^  can  now  be 

itodifieo  as  follows  to  obtain  corrected  values  that' are  closer  to  D 
and  Yj  for  =  Oj  ^ 


where 


f,iaod 


r  -| 

,  expr  D  ^ 

f , expr  J 


{VII-19) 


^  -  D 

^  f « 0  f,  expr 


{yiI-20) 


“f.expr  “"P'-'Au’ 


(VII-21) 


where 


v7" 


a  -  V 

^  f»X=0  -f,  expr 


(VII-22) 


(VII-23) 


expr  ®2 


{VII-24) 


^f,sod  ^f,!aod  ^  used  in  Eqs.  {VII- 17)  and  (VII- 18) 

to  obtain  new  afproxiisations  for  f  and  hence  The  whole  procedure 

can  be  carried  tiirougfa  as  often  as  needed  to  cotain  the  desi’-ed  accuracy. 

^f,s>d  \sod  rapidly  to  axid  Por 

exar^le,  in  the  calcalaticsis  for  Core  510,  only  two  iteration  cycles 

were  used  at  cost  by  the  cceauter  progran;  the  t.Hird  iteratioTa  cycle 

produced  a  flange  in  o  of  less  than  O.i  percent  of  D.  .  The  total 

i.,  expz 


10? 


h 

c  p 

A  C!  CVZJ— 2o) 


i03 


''•'“'•'Ml*™, 


Friction  Factor  Calculations 


The  relation  used  for  the  calculation  of  the  friction  factor  is 
fran  [14]* 


A  V-  r 
_  c  1 

A  V 
m  L 


ZiP  X  2g 
c 


=^v. 


-( 


K^tl-p  )- 


{VII-23) 


where  the  subscript  1  indicates 


* 

X 


0 


and  2  indicates 


* 

X  =  1 


and 


V  =  (V  -r  y  )/2  (VII-30) 

n  1.  ^ 

^e  tensperature  used  in  the  calculation  of  is  the  tine-isean  teapera- 

* 

ture  of  the  fluid  leaving  the  heat  exchanger,  T  (x  =  1) . 

m 

Periodic  Technigufi-Coaputer  Prograa  for  Data  Reduction 

Olie  data  reductian  progran  used  in  tiie  periodic  technique  to  cal¬ 
culate  the  heat  transfer  and  flew  friction  characteristics  is  described 
in  this  section,  Ihe  pixgran  is  written  in  the  Fortran  IV  language  and 
was  used  on  an  ISM  360/67  co;3j«ster.  It  is  bajxid  on  [4}  but  has  been 
extensively  refined,  ihe  prograa  has  been  successfully  used  for  a  test 
range  of  0.4  <  <  36.6.  It  is  designed  to  be  used  where  one  ter^sera- 

ture  record  is  taken  upstreaa  of.  the  heat  exchanger  test  core  and  the 
other  is  taken  downstreas  of  the  core.  The  prograa  is  is>t  sensitive  to 
the  point  in  real  tine  used  as  a  starting  point,  but  the  rastreas  and 
downstreaa  records  nast  occur  sinultaneously. 


Input  romata  and  >kaaenclature 

Input  is  in  the  fom  of  data  sets.  A  data  set  consists  of  test 
core  paraneters  {Cards  i  through  3)  followed  by  individual  run  infoma- 
tioa  CCarv*^  4  thnsiOT  8)  ail  using  the  sase  set  of  core  paraasters. 

Any  msrber  of  data  sets  nay  be  run  in  series. 


IGS 


I 

I 


C^Xd.  Contains  fchA  rw^av-alTi  _ _  «  ^.  .  . 

- -  -i-wucixxcatxon  xn  columns 

1  through  7o.  h  nonzero  digit  in  column  80  is  also  required. 
This  card  serves  a  dual  purpose,  it  initiates  a  new  data  set, 

and  it  causes  a  svnnmary  of  the  previous  data  set  results  to  bl 
listed  should  one  exist. 


Card  2 


Card  3 


Card  4 


This  card  contains  the  test  core  parameters  and  geometry  infor¬ 
mation  for  the  series  of  test  runs  that  foJlow.  Each  entry 
should  contain  a  decimal  point,  ihe  columns  allotted,  variable 
names,  definitions,  and  dimensions  are  as  follows; 

Col.  1-8  ALP  -  core  a^ha,  ftVft^ 

9-16  POR  -  core  porosity 

17-24  Afr  -  core  frontal  curea,  in^ 

25-32  I£KdH  -  core  flow  length,  in 

33-40  -  core  weight,  lisa 

41-48  SHTP  -  core  specific  heat,  Btu/dtaa  ®f) 

KC  -  core  K 
c 

57-64  KS  —  core  K 
c 

65-72  COND  -  core  coTsductivity,  3tu/{hr  ft  “F) 

This  card  contains  the  dri*  bulb  tei^ature  ana  the  wet  bulb 
temperature  for  use  in  calculating  the  airstream  humidity. 

Bach  entry  should  contain  a  oeciaal  point. 

Col.  1-e  TiB  -  dry  bulb  tusserature,  “F 
»-16  IW3  -  wet  bulb  teiserature,  “F 

This  card  contains  t.ie  individual  test  run  description  in 
coiusns  i  throa^  76.  Colusi  80  is  left  blanjt. 


Card  3 


..ol,  1—2  K  -  nc^jer  of  t^perature  points  loich  are  read 
for  each  terperatare  record  is  tiJis  test,  so 
decimal  point  is  needed,  bst  right  adjust. 


no 


9-16  Rj  -  reference  jujKiti*::  ievei  is  riiiiivoits  icr  i 

copper-constant<in  tJierBs>co?;sle  vith  an  ice  bafe. 

17—24  ?ERIS>  -  the  period  S  of  the  individual  Zi&sz  mn, 

o 

sec. 


Card  S 

Col.  1-30  T(I)  -  ten  of  the  vipstreas  temperature  data  points 
which  are  read  off  the  upstrsas  record.  ?.s  nany  of 
these  cards  are  used  as  is  repaired.  coissss 

core  allowed  for  each  tsrseratare  reading  in  nicro- 
voits,  with  ho  decinai  point  required,  •srt  each 
entry  nust  be  right  justified  in  its  3  c^iunns. 


Card  7 


•This  card  is  sinilar  to  card  6,  except  asanstrear? 
readings  in  sicrovoits  are  entered. 


Card  8  ahis  card  contains  fluid  flow  rate  infomation  for  each  insi- 
"Tidual  test  run.  Decisal  points  are  repaired  in  all  entries 
except  ICSiF, 


Col.  1 


ICSJF  -  nusber  of  orifice  plate  a^d.  OrSnics  plats 
nsFber  1  is  1.301  iiKhes  in  dianeter  and  nanher  2 


is  2.201  iixhes  in  diameter. 

2-8  E£^  -  crifine  pressure  drop,  in.  water 

9-16  iHO?  -  core  pressure  drop,  is.  water 
17-24  3jS  -  anbient  pressure,  bareneter,  in.  sercury 
25-32  ?C  -  deprsssicHi  of  core  islet  pressure  below  nnbient. 


m.  water 

32—40  ?0  —  depression  of  orifice  inlet  pressure  belrst 

At.,  w,  m.  water 


To  corpute  sore  than  c»e  test  run  with  the  sasas  core.  Cards  4  throng  S 
are  repeated  as  often  as  required.  If  the  core  is  ch-c^d.  Cards  I 
throo^  3  2s?st  again  be  used.  At  the  end  of  all  r^ms,  a  slash  card  wirh 


diuit  in  colum  SO  is  usad  to  cai^e  a  sunnary  of 


a  ocnzerc 


A  Fortran  Sourca  Listirjg  c*’  the  Perio«lic  -Transient  Technique  Datra 
Reduction  Progr^  is  presented  in  Table  VIX-3* 

Additional  Informaticn  on  the  Use  of  the  Fericdic-Test-Technique 


Data  Raduction  Trograat 

The  prograa:  consists  of  three  parts;  tiie  nain  prograa,  the  sub¬ 
routine  COKVSS,  and  the  siibrou  -  ine  CONDCT-  The  Kaia  prograra  does  the 
bulk  of  tJie  calculatione  sad  all  '-f '  the  input/output.  Stsbroutii^  CONVES 
cpftverts  thersc»cotfflIe  voltage  to  degrees  Fahrenheit.  Subroutine  CONCCT 
caicuiatea  rising  inputs  or  ^  and 

A  constant^  CHART  iiP  included  in  the  prograr.  as  a  scale  factor* 

It  provider  flexibility  and  convenience  in  the  type  of  recorder  used  in 
the  thSrJBocouple  voltage  input  Cards  C  and  7  and  the  number  of  thernso- 
couplee  used,  ultimately  the  scale  factor  is  used: 


‘‘actual  'input 


X  CHART 


Physical  properties  at  the  solid  core  material  and  the  physical  properties 
of  air  used  in  the  periodic  data  reduction  are  identical  to  those  dss- 
cribed  Jn  this  appendix  -under  the  Single-Blow  Technique. 


Description  of  Subroutine  CQNDCT 

The  subroutire  CONDCT,  Table  yil-3,  calculates  and  for  values 
of  i  and  and  is  used  in  the  solution  for  the  periodic  technique 

witii  finite  wail  longitudinal  conduction,  Eqs.  (V-25)  and  {V-26) .  It 
is  entered  by  means  of  a  fortran  CALL  statement,  A  description  of  the 
parameters  in  the  SUBROUTINE  statement,  follows; 

LNTU  -  product  XN^^ 

p  -  nondimensionai  period  ^ 

CEP  -  e,  for  X  =  0  or  1/(1 

f  2 

OET  -  62  for  X  =  0  05^  +  ?  ) 

EO  “  used  in  Eqs.  (V-25)  and  (V-26) 

ET  -  e..  used  in  Eqs.  (V-25)  and  {V-26) 


Table  VXI— 3  p35®  I 


C 

C  ?C«!031C  TRtHSIENT  Tgsi  rECMf<IOi£  CIT*  ^FfclCTIOH 

C 

1  itEM.  lsENGrHt»7Uf;iTUfttXCk<E(MU 

2  OtHENSrosc  ALPmi^l  t3ET(l9l. 

3  OtHENSION  BFt lSOI*AFt ISO! tBI 11.21 tA(li»21 }*T llHOI 

4  OtHENStON  £Vi20}  ttC(2CltN<i:201«S>E{20(«PSAt20)«>.KTAt23i7A*.Ct2ClfAST 

t>K  ( 201 .  SP SSI  201 .  Si^l  A{  20)  .  SRC«20)  .  SSTPKS120}  *UFI2C  1  * AiAHt2C  3 

5  ?R«>q.710 
TCC^sO.C 

7  CM*RT-'0.5 

8  0C=0 

S  BAR-0.0 

C 

C  READING  IHE  TEST  CORE  PARAJiETERS 

C 

10  1  REACi5*49l  At.PHtINACT 

11  F0«K4T<1«)A4.3X,1U 

12  ifiinact.eq.oi  go  to  3 

13  IFCNC.NE.di  GO  TO  43 

S-i  4*  DG  47  JO- It  19 

15  47  getlJO)*ALP«)JC) 

16  REA!>(5f483«EI!0>60IRtF«POR*AEtt«LENS?H(HTtSPHT«KCsS;E  fCONO 

If  A«All,»AfR»U.0-Pt«r/l44. 

Ig;  AC»FO«*AFR 

19  REA&{<-.480fTDB.TJa 

20  CW«WT*5PKT 

21  t7»At?*AFR9tEH&THyif28. 

22  STP=IAC/ftr  l*PR’**l2./3.)/l44. 

23 

24  CO  TO  I 

C 

C  2=401  NG  THE  DATA  fOR  EACH  INOIVIO'JAL  TEST  CASE 

25  3  REA015.491}N»RJ»PERI00 

26  NCiNC  *  I 

27  00.2  i-.l,2 

28l  !iEA0I5,580MT<I»ii-i.»« 

29  GO  TO  1 41,  421  ♦! 

30:  41  MRirS(6i4881 

?l  WRITEt6,54) ,AtPM 

32  54  F0R)*ATaK,19A4/l 

53  X?rTEI6,48tl 

34  GO  TO  26 

35  42  i{RnEt^,489) 

56  26  CONTINUE 

37  KRITEffc»482»IJtT:J!>J»l,Nl 

33  ft-*r*R4 

c  ,  - 

C  CONVERTING  HltLlVOtrS  TO  DEGREES  ? 

C 

39  00  19  1*1, N 

40  Til )«T *11 •CHART* 9) 

41  19  CAtl  CCNVRS(T(n,TCaN) 

42  CALL  CCNVRSfRJTtTCONl 

43  HRITEI6, 432)1 j,Tt45 ,J»l,N) 

44  WBI7Et6,499) 

45  OSlG*x./FLOAT{N-l) 

46  PI*3. 1415920 
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C 

c 

c 

*7 

*8 

A9 

50 

51 

C 

c 

c 

52 

53 
5* 

55 

56 
5T 

58  30 

55 
60 
61 
62 
63 
66 

65  40 

66 

67  21 

68 

69  27 

70 

71 

72 

73  2 

74 

75 

C 

C 

c 

76 

77 

78 

79 
-80 

81 

82 

83 

C 

c 

c 

84 

85 

C 

C 

c 

86 

87 

88 

89 

90 


?I60!KC  THE  E0U8IER  C&£»=FIC!ENTS 

00  21 

sroo.o 

Ni«HM*l 

N2*H-l 

eR*2  1G*FL0*T  tNH 

TRAPeZOIOAL  R^E  fNTEG8ATI09 

4F<  llx^^rt  11  *S!K<8*l/2.0 
8F11 >«Tti 1 •casieRi/2.0 
00  3C  J«2«f< 

SIG>SIG40SIG 

8Ri? -0*PI*SlC«FL06r ^Mll 

AFIJ}*T(JI*SIN(ait} 

3F  K*  }«Tf  J  J  {  88 1 

AFIN}«AF(NI/2. 

aF{Nl>BF(H)/2. 

BfKHtiLt-0.0 

*(NN*LlsO.O 

00  AO  K^ltH 

e<NtlfLl-aiHN.Ll*S=(Xi 

A{m«6l*AlNNtL}«'AF{XI 

BIM«.Lt-8fNN>(.l/{0.5*FL06Ttli2}) 

A<fW,!,»*A{*{»l»tl/f0.5*Ft0*TI«(2ll 

00  27  KI*t,a5 

F8tKn«Kl>l 

<lX!rE{6t495Hl(Mf  11  *A{  !,LI.I«t«ll} 

HBlTEi6.496}{MM(n«S{l«(.).l-i«lll 

WtlTE(6*510|l 

CONTtHUE 

r*«AN-B(l*2}/2. 

HU-0.C4254 1 5.  lE-051  •<  TH£Ali-40.3l ♦l»E-07* <  THEAII*40,3I**2 

CALCULAHMS  THE  FLOk  FA8AMETEIIS 

llEA0f5«568)IORF»HEA0r0«06t0AB*FCt8O 

WR I TE  <  6 1 5091 8 AR  «81T  ^TCOM.TiCANt TOBf  TMB t ALF  «  KC* KEvHT 

Pl«8AR*70.73-PC*5.204 

P2»Fl-08(:^*5.204 

F3»Pl-P0*5.2O4 

P4*P3-HEA0*5.204 

X»IP3-P4I/P3 

Tl*459.67»^Tn£Ak 

DENS  ITT  CORRECTICSJ  FACTORS 

H»0.0196919-0.030679526*T1«*  8.916A7E-96*TM86«2-9.000232* 

lfrCB>TMB) 

XV«(l.0^l,607*HI/ll.0*HI 

SPECIFIC  HEAT  OF  AIR 

CP»l0.87*H»t.0001*0.24 

Vl»fXV*53.?4*Tll/Pl 

V2»Vl»Pl/P2 

V3»Vl*Pl/P3 

V4*V1»PI/P4 
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c 

91 

92  ll 

93 
9* 

95 

96 

9?  Zt 

98 

99 
100 

101  5& 

102 
103 
<04 
IC5 
106 

107 

108  560 

109 
lie 
111 
112 
123 
119 
115 


SSl.£CTI%  T}^  FlOU  COEFFICIEHIS 

60  lO  111.22), lOKF 

01-40.6053 

9£TA>0.62S2 

&{ 60-1.301 

OC-0.012T 

00  TO  50 

CF>Q.110Q 

8ETa*0.7193 

o;:  -  ••  1 

^  <4I*0.35*3ETA**«-)*x/i.4 

l*4*359.*C?»  t0I*0*92jeS<^  tHEAD/Vi  fr 
90-48.  9Wi/<FI*0160*m)i 
CF-tCF^ik:*!!  .OE«C4/KOI«aO.S 
Wfs?58i«CFr{£')i6o**2}*SQilTIHEAO/V3}*r 
WS-yH/3600. 

L*IB06-CO!i09AWAt.L*12-/tL£MGT)^VH«CF) 
F0»I6TI  X  I,  F?.2,  SF8.2) 
i)Hs4.*AC«LEM&rK/ 147*1728. 1 
M6FVfKI«eH*0H*l  44.c/()tU*ACS 
6ftITEf6,SlC) 

K*tTE(6.50S) 

WKITEI 6,506) 

:^I7FI6,507| 

6*6K*i44.0/6C 


C 

c 


116 

117 

118 

119 

120 
121 
U2 

123 

124 

C 

C 

V* 

c 

c 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134  71 

135  72 

136  73 


CAlCOlAItKi  DC  FAftriOn  FACrOX 

V***IVl«¥2»/2.0 

05-65*144. 0/AC 

^”^^380j**5.204*2»69-32.2)/16S**2*VH—C  mC*1.0~FOA**2l»>  fi*tw9/wt  i 

J  ******^  •*^®^*^****®*?  ¥*.»  *tt  AC,  AF*  .IgKCTMfCW,  Oh,  FWIOOr 

*«IT£  16,5101 
WKtrEf6,699) 

M1TEI6,494} 

WKITEf 6,499) 

MRITEfe,698) 

USIMG  THE  roURlEX  COEFFICIENTS 

PHASE  SHIFT  CCMPMlSON 

00  24  J J-1 ,5 
XFfN.E0.3l  60  70  35 

M^-FLOATIJ  J)  *12  .*8 1»C||  }/f  PEA  I00*MS«CP  » 

ttNN-A  rAN21 1 A1  JK,2 1 /tC4K,2)-Al  JA  ,  ll /8C  3A,l| )  ,  11  .♦M4K. 
12)*AlJK,l|/f8f4K,ll*01JK,2)))) 
tFffERWI.&T.O.Ol.AW.ftfIS.Ot.S.OO))  CO  TO  72 
XFtfEAHN.LT.0.0).AN0.ltlNS.6T.4.25))  SO  TO  71 
IFIEANM.LT.Q.O)  GO  TO  72  •  ou  n 

(»  TO  73 
EAPK-EANNFFI 
ERHN-EANN*FI 
CONTINUE 
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S3? 

IBS 

!3$ 

ISO  ?4 
14S 
142 
1« 

W 

145 

145 

14? 

145 

145 

150 

151 

152 

153 

154  ?4 
C 

c 

155 

154 
157 

155  ?t 

159 

160 
161 
162 

163 

164  - 

165 

166 
16? 

16S 

169 

170 

171 

172 

173 

174 

liS  77 

176 

177 

178  88 

179 

180 
181 

182  89 

183 

184 

C 

C 

c 

185 

186 

187  35 

138 

189 

190 


lfSf«0 

056*0,9 

a>5F^0 

C£4II*C«I9I*0C? 

V£*iWi0*S/€®Ul 

lM7£8l|if<.i.T.fl.O}}  SO  TO  8 
l./iQSTlTSUdKl  ,0  } 

RCC*9E8tC»/ IFLCATI 
!»TU*F£€*st!iS 
IFiJJ.KC.l}  TO  76 
Lfe?UH.8MB9A«IITU 

C*U  COS>CT(Li(ru«F££,060*(%T,£0»£TI 
046«I0£T"^TI«KTU 
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fSje  sastiod  ficl  is  csed  to  fioi  tfee  issts  o£  r?-2?j 

asd  |?-22).  Inis  is  ass  iterative  ?!iet.!»d.  jaitial  ^s^ssss  ^seS  ais  for 
2  2 

c,  =  i/{l  +  t  3  ®,  -  */0-  -9-  r  3,  tae  i  =  O  ca^.  Ccnveraes^ 

is  very  xz^id  ard  iteraticti  is  st^^ed  vtea  ^:e  c^sectiea  besanes  less 
tisaa  0,1  percent  of  or  €,-  'Ssis  ©scarred  cariag  trs  third  itaratics 
is  tlse  excregg  ca^  with  is.  =  3.5  aad  f  =  S. 

Periodic  ^'echaigae  —  Outrxit.  of  Test  basalts 

?uo  foszs  of  the  oataat  are  "crided  by  *he  aerioaic  cectgiista  data 
sessTitics  progrzE.  tor  cacis  ra:^  s  basic  listiog  of  thesraccssssslc  voltaye 
the  correspordirg  teJi^raUireSj  the  POarier  Coefficiests,  core 
pas^etsxs,  fio»  rate  paraseters,  alcrg  vith  beat  traasfer  saso  f lose 
friction  zesalta  sst  provided.  At  the  ecd  of  erach  data  set,  a  sasnsary 
of  ail  the  xaas  for  Sie  is  listed.  Ihe  ^osecdatcre  for  this  tistlrg 
is  ^vcs  is  5able  yiI-4  aoS  tbs  zssslts  of  testiso  c>sr»  SIO  are  presssted 
is  Table  VH-S. 

Table  Vn-4 

5K^SE2«G^'»SHl  .mSlCSI£-  TSEHSI'^  HSSSI^ 
w  rloaerate* 

G  Ssdss  velocity,  liat'ibr  ft^) 

la  i^eysolds  nsssscr,  diaess  jg^lsss 

^SiSCS  stricd  of  osciliatits  8^»  sec^ 

o 

IfiSSDn.  X,  heat  traasSir  ccacscticn  paraaetsr 

PSI  f,  ncsdiEsi3sicciai  period  of  os^dilation 

sro  H.  ,  dis!?ssioaiess 

SC.  ^  tiee  constant  of  the  core  solid  naterial,  see© 
j  Colbom  3  factor  »  diKSsiosless 

5  Friction  factor,  f,  disehsionless 

Note:  The  parameters  »I,  RTU,  iCSJ,  and  J  are  presssted  twice. 

5Sie  first  set  for  the  asplitcde  attenuation  sethod  uses 
£<J.  {vii-17;  to  co^tB  f.  The  secijnd  set  for  the  ^ase 
^ift  uses  Eq,  (Vil-iS)  to  corpute  y.  Althcsugh  bod;  are 
presented  for  eaadi  run,  oiily  the  results  of  tne  recos»naed 
sethod.  Fig.  7.3,  should  be  censidered  useful. 


119 


PKRIOntC  rCCHNIttilt*  f*xr*F(UMPNTAL  KKSUI.VS 


T 


^  o*  y  — ^  ^  ^  s 

ii. 

.  SOO 

Cw^cc.  SiC-C  wc.ee 

—  ^ 

C?  ^  y 

^af  C£ 

S  O  C  3^  ^  ^ 

«•  -C  ?-  C.  £  c 

-5 

?!r«s,-«  — «»-c  c  eocc 
C  CCCrC  C'CCCCCC 

»  w  ©  ■« 

^  ■©■  •  •  *• 

Cc  OCCcCCOCwC 

•  »  *  O'  .  • 

c  oo?^  ewc: 

a* 

£qr  <C5 

^3  0 
s.  war 

M  ©  K 

.c^  «— C‘Crcc£«3r£^ 

^  ^  m  9>  X 

Ci, 

<  »• 

«A 

cr  ^ 

*/!*- 

<  V?  C£ 

■«  z 

^"?  0“  1^  SSi  -aff 

mm 

mm9^  ^  ml  ^ 

m 

rjeu 

>-  M  sz 

s\:  »>i 

o 

<2  O  w 

KT- 

^<9^^  OSC-Cf-SP^,^  — 

era,  ac>  £ 

m. 

ae  cc:  ac  -^Sif 

&.  &  o  £: 

•—c  ?-  -«£©  'OOS-oate 

j  ■< 

« 

mk. 

o^cfcoc  ooewoc 

mm 

C  O  O  Cr  O  c- o  O  o  O  O  O 

■<  • 

j*“  rw  ^  ^  0» 

9* 

->  3;  W 

ass  £££S£35££r»P- 

mtSi 

Z  wss 

o 

& 

V  O  X  w  ^  ^  O  w  O  O'  O 

»<i 

1 

t 

^  S£ 

■c  , 

^Ot».^P-£.'Z~©' 

u.  t  i 

is 

kO-  -C  £&’:<&  0^£ 

o 

ejf-S£  a 

z 

'«‘-J  'Z 

-y 

£.Pi  £iP4_c^P*Cf'*-'««^-Ss, 

IM 

-  »-’V  ^ 

X 

^ 

z: 

.-tt.3  3 

ZCf»-  »- 

•j 

a.  a: 

j*ii-sr.^®f^-««6sP£mew 

w 

e^oc-cco 

Z9l 

«  «  »  o  o  -«  o 

<-*^-  --f«  aOM 

•<' 

^  ^  -«■  rt  «  £•<•«•««»•—»« 

o 

C  •  9^  •  •  ■# 

O  <M  '4  O  «  • 

£ 

5P«P  PSOtt-'C^*-*  «^**»««fV 

*  ^  o 

fte-“'-«»-o  coooooo 

o 

S 

O  O  O  O  C- O  C-C  O  £  C  o 

z 

oc  £Ooc-coceo:> 

« 

«  ?  a  c  K  ■  « 

•« 

c 

OOOCCC'CCOOCO 

z 

<  c. 

«  z 

a: 

C  X  ttT  ~ 

w  « 

C:0£'9 

c 

»-»-*-  • 

—  c? 

•CC  ^^-Kir^O<^—rv^-e 

U>  X  o 

ac  w 

•~r~  £  z 

u» 

.©ePCiiC  P-'-C© 

w  w  a.  _j  c 

S  J  ««  3  V? 

■«  >«  £ 

rrs  p*  Pi  PS  ^  — 

£  3  X  U  u. 

COa.ac»-  cc 

oc 

>  —  jC.u;cjg 

Z 

^  P-  --  .c  £  tv  £-«!-•>©--» 

X  B-  <  U3  3C  u 

— <  -C^-  tw  tN4  PS  ^ 

- 

O 

©  ev  p-  tv:  *V  >C-5P  £r»»-  O  O'  .<• 

©©•«r. OC  -*s:sv«>»w»v:o 

•v. 

^  ©  S'  C  ®  ®  C-**  'p  © 

3 

»«  —  P-Pif  '©'©©  ©•*:«'« 

K 

s:»-£®<.s's»--ca  — p-o 

-*svO»^Oacfv^—  »t*'s a  ^- 
£  O"  £  T-  ?- 


120 


D?elassifi^ 


DD  ,^..1473 


Compact  heat  exchangers,  testing  techniques 

Periodic  technique  for  testing  compact  heat 
exchangers 

Compact  heat  exchangers,  design  data 


